











Use of Callendar's “Radio-Balance” for the Measurement 
of the Energy Emission from Radioactive Sources 


W. B. Mann 
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1. Introduction 


In 1910, in a paper that seems to have caused 
little interest in the intervening vears, H. L. Callen- 
lar [1], deseribed a calorimeter in which the quantity 
ff heat to be measured was to be balanced by the 
cooling from a Peltier couple. With this calorim- 
eter, which he called a radio-balance ° 
gave a null reading, he made measurements of the 


because it 


solar radiation, and he also carried out a somewhat 
rudimentary measurement of the energy emitted in 
unit time from a sealed source of radon emanation 
lent to him for the by Lord Rutherford 
The results obtained with the emanation, although 
qualitatively interesting, were not of value 
from a quantitative point of view 
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The radio-balances are shown, from left to right, in chronological order 
radio-balance is that with four Peltier junctions 
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later Hoare [2, 3] used a radio- 
balance to determine the Stefan-Boltzmann = radi- 
ation constant. This appears to be the only other 
recorded application of the radio-balance, and no 
other reference to it in the scientific literature of 
the last 23 vears has been found. It seems, how- 
ever, that this method of determining the energy 
dissipation has much to commend itself, and experi- 
ments have recently been carried out to see whether 
or not it might be feasible to apply it to determine 
the rate of energy production by radioactive sources 
For this purpose, several radio-balances were con- 
structed. The outstanding merit of this form of 
microcalorimeter is its extreme compactness, the 


Some 20 vears 


first simple version used in the experiments here 


described being contained in a cylindrical copper 
box only 2\ in. high and 2 in. in diameter, whereas 
the final version was contained in what 
tially a evlindrical copper block 24 in. in height and 
fin. in diameter. These are shown in figure 1 
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2. Theory 


In the racdio-balance the energy incident upon, or 
transmitted to, a small disk or cup is balanced by 
the cooling effected by a Peltier junction attached 
to the same disk or cup. Callendar used both a 
disk and a cup radio-balance as a pyrheliometer, 
and Hoare used a cup balance to measure the Stefan- 
Boltzmann radiation constant 

In the present experiments a racdio- 
balance was used, consisting, in its final version, of 
two copper cups '%¢ in. long and having internal 
and external diameters of 0.099 and 0.130 in., re- 
spectively. Each cup was mounted on a six- 
junction constantan-chromel thermopile by means 
of an equatorial ring, as illustrated in figure 2. The 
lower junctions of each thermopile dip into silicone 
vacuum oil contained in two small cavities in the 
lower copper block constituting the container of the 
radio-balance. A constantan-chromel Peltier couple 
is soldered to the base of each cup 

The radio-balance may be used in two ways 
Each cup may be used singly, or they may be used 
together in opposition. In the present experiments 
the interest has lain in the measurement of the 
energy emission from radioactive sources and not 
of energy flux from an external source of radiation 
Such a radioactive source, as, for example, a radium 
needle, may be placed in either cup A or cup B of 
the radio-balance and a current passed through the 
Peltier junction in such a direction and of such 
magnitude that the cooling due to the]Peltier june- 
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Elevation drawing through the vertical azis of sym- 
metry of the third and final radio-balance 
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heating due to the 


tion shall exactly balance the 
radioactive preparation. The criterion for sueh 
balance is obtained by reducing to zero the deflection 
of a sensitive galvanometer connected to the thermo 
pile For zero deflection, the temperature of the 
cup is restored to the same temperature as the bas: 


| of the thermopile, and there should be an exact bal 





ance between the heat energy supplied to the eu; 
by the radioactive source and the heat removed by 
the Peltier cooling. In this condition the energy 
I], supplied by the radioactive source is given by 


H=PC—C'R, 


where 7’ is the Peltier coefficient, C is the current 
passing through the Peltier junction, and F& is its 
resistance 

In order to determine the resistance of the Peltie 
junction, a subsidiary experiment may be carried 
out, in which a dummy source is substituted for th 
radioactive source, in order to simulate external 
heat-loss conditions and the thermal conductivity 
along the length of the cup, and the current is again 
adjusted to give zero deflection of the thermopik 
galvanometer. In this condition 

0—PC—C*R, 

or 


R P/C (2 


By differentiation of eq (1) it is clear that the max- 
imum cooling effect that can be obtained is 


H=P*/4Rn, (3 


In the second method of using the radio-balance 
the thermopiles of cup A and cup B are connected in 
series with one another and with the galvanometer 
so as to act in opposition to each other. Thus an 
equal rise in temperature of both cup A and cup B 
relative to the thermopile bases will give a zero 
galvanometer deflection. The Peltier junctions are 
also connected in opposition and in series, so that 
the same current flows through each but will cause 
a heating in the one cup and a cooling in the other 

In all the radio-balances that have been constructed 
in this investigation a small residual galvanometer 
deflection was almost invariably detected, indicating 
a slight temperature difference between cup A or 
cup B and the thermopile bases, with no source in 
either cup and no Peltier current flowing. This 
could be due to a very slight difference in tempera- 
ture between the top and bottom of the copper 
block container. With the piles of cup A and cup 
B connected in opposition, however, the “zero de- 
flection” was invariably negligible, any slight rise in 
temperature in cup A being compensated for by a 
similar rise in cup b. 

In this second method of using the radio-balance, 
which was also proposed by Callendar, another 
advantage derives from the fact that the Joule 
heating of the Peltier junctions is also eliminated. 
To describe this method, assume that a radioactive 
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source dissipating energy HW" watts ts placed in cup 
\ and that another dissipating W’’’ watts is placed 
in cup B. If a dummy ts used in cup B, then W”’ 
is zero, but this is not necessarily so, and the method 
can im fact be used to compare two nearly equal 
sources Without involving the use of Peltier cooling 
or heating at all. This appears to be a new and 
very accurate application of the radio-balance. as 
will be demonstrated by the experimental results 
obtained. This application would certainly have 
been most apparent to Callendar from his derivation 
of the relevant equations; the need for such an appli 
eation did not, however, exist to the same extent 
more than 40 vears ago 

it will be assumed that W’ is greater than W’’ and 
that a current Cis passed through the Peltier june 
tions in such a direction as to give a cooling of cup 
( and a heating of cup B. The value of ( is then 
adjusted to give an approximately zero balance of 
the thermopile galvanometer. Suppose, however, 
that the residual galvanometer deflection is d The 
sources HW’ and W’’ are now interchanged between 
cup A and cup B and the current C is reversed so 
that cup A is now heated and cup B is cooled. As 
slight differences will exist between cup A and cup 


B. a new residual galvanometer deflection, d., will 
bye observed Under these conditions Callendar 
showed that 

WwW’ —W”’=2P¢ d,—d,)/s, (4) 


where s is the galvanometer-scale sensitivity in mil- 

limeters per microwatt transferred from cup A to | 
cup B, and P is the average of the values of 7? and | 
P, for the individual cups. In the final version of 

the radio-balance used in these experiments, s was of 
the order of 0.7 mm uw when the current through the | 
galvanometer was reversed. That a reversal of the | 
current does, in fact, correspond to a transfer of | 
energy equal to 27’C' from the one cup to the other | 
is clear from the energy-dissipation-level diagram in 
figure 3, which refers to either cup A or cup B. If | 
cup A gains 2PC, then cup B correspondingly loses 

















2PC. The value of # in figure 3 would be that 
appropriate to cup A or cup B, 
| 
Pc+ ctr | 
-——— — ee ee ee ee Pc 
2PCc 
° 
- pc+c®r 
—_— ee ee eee eee ee ee ee - Pc 
Ficure 3, Energy-dissipatron-level diagram illustrating the 
elimination of Joule heating on reversal of the current through 
the Peltier junction of either cup. 
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3. Experimental Procedure 


The cups of the third, and final, radio-balance were 
constructed from copper and were 0.099-in. internal 
diameter, 0.130-in. external diameter, and '%q in 
long. The cups were supported on the upper end 
of a six-junction chromel-constantan thermopile by 
means of an equatorial annular ring with six small 
holes. The junetions of the thermopile were in 
sulated from each of the cups by means of six disks 
of rubber hydrochloride, 0.0002 in. thick, which 
were cut out by means of a paper punch. The 
junctions were then sealed in position by means of 
a small amount of glyptal. 

The Peltier junctions are secured to the lowe 
ends of the radio-balance cups by means of soft 
solder. The Peltier-junction assembly consisted of 
a framework of four wires in the form of a Cleopatra 
needle silver-soldered at the apex. A constantan 
wire and a chromel wire, each 0.020 in. in diameter, 
along opposite corners of the needle formed the actual 
Peltier junction. The other two wires formed a 
thermocouple, by means of which the lowering in 
temperature of the Peltier yunetion could beymeas 
ured. These consisted of constantan and chrome! 
0.010 in. and 0.012 in. in diameter, respectively 
When assembled each cup stood upon the apex. of 
the Cleopatra needle in the manner indicated in 
figure 4, 

The thermopiles were constructed with’ the help 
of a number of lucite annuli of %-in. external diameter, 
se-in. internal diameter, and % in. thick. ‘Twelve 
small holesjwere drilled through each of these annuli 








Figure 4. “‘Cleopatra-needle’”’ assemblyjof Peltier junction and 
thermocouple for measuring the temperature of the Peltier 
junction. 











parallel to the axis and at 30-deg intervals.  Alter- 
nate O.010-in. diameter constantan and = 0.012-in 
diameter chromel wires, 6 lengths of each, were 


now threaded through the 12 holes of 6 of these 
annuli in series. The annuli were separated by some 
2% in. from each other along the wires and secured 
to the wires by painting with glyptal cement. When 
the glyptal hardened the wires were cut midway 
between each anulus, and the basic units for six 
thermopiles were obtained. A thermopile was then 
simply obtained by hard-soldering adjacent constan- 
tan and chrome! wires at each end in such a manner 
as to give one continuous circuit 

The cup and thermopile assemblies were finally 
mounted on two lucite disks in a copper disk mount 
in the manner illustrated in figures 1 and 2. 

An intermediate type of radio-balance, illustrated 
in figure 1, was also constructed, in which four 
Peltier junctions were bound with nylon thread to 
the cup, 2 above and 2 below the thermopile, to 
investigate the effect of distributing the cooling 
effect However, apparently due to the fact that 
no part of the cup is ever very different in tempera- 
ture from the surroundings at balance, this type 
had no advantages over that of simple construction 
and was, in fact, less satisfactory in use. 

The Peltier current was derived from a 6-v battery 
in series with a 1,000-ohm variable resistance and a 
standard l-ohm and standard 10-ohm coil. The 
value of the current flowing through the Peltier 
junctions was then determined by measuring the 
potential drop across either of the standard resist- 
ances by means of a potentiometer. The current 
from the thermopiles was measured by means of a 
high-sensitivity galvanometer of nominal sensitivity 
equal to 0.000029 yaimm. The thermopiles were 
connected either individually to the galvanometer 
or in Opposition to each other through the galva- 
nometer by means of mercury switches. A reversing 
switch made from three-contact mercury switches 
enabled the current through the galvanometer to be 
reversed. To eliminate extraneous thermal effects 
copper wire was used throughout, and the connec- 
tions to the mercury switches and the switches them- 
selves were immersed in an oil bath. An 18° to 28° 
C thermometer was used to measure the temperature 
of the copper block containing the radio-balance. 
This enables the variation in P with temperature to 
be calculated, using the thermodynamical relation 
that 

P= TdE/dT, (5) 


where 7 is the absolute temperature, and £ is the 
thermoelectromotive force. 

In the present series of experiments the perform- 
ance of the radio-balance has been very extensively 
investigated by using ordinary commercial radium 
preparations contained in platinum-iridium needles 
having a wall thickness of 0.5 mm. 

The value of P may be determined in a subsidiary 
experiment by measuring the variation of thermo- 
electromotive force with temperature of a thermo- 
couple made from the same wire as that used to 





This in fact w: 
done, the thermoelectromotive force being measure 


construct the Peltier junctions. 


at 25° and 50° C. In practice, however, a slig! 

temperature gradient of a few hundredths of a degr: 

will exist along the cup between its base, to which th 
Peltier junction is attached, and the equatorial rin 

of the thermocouple junctions. This means tha 
the cooled cup will gain a small quantity of heat fro: 

the surroundings, and the heated cup will lose a co; 
responding amount. The net effect will be to cd 

crease slightly the effective Peltier coefficient | 
is better therefore to measure the Peltier coefficien: 
under the conditions of the experiment in the radio 
balance itself. 

For the purpose of determining the Peltier coeffi 
cient in this way, two pairs of small resistance coils 
were wound on platinum wire to simulate the thermal! 
conductivity of the platinum-iridium needles. Each 
pair of coils consisted of a main resistance (2?) with 
consisting of only a 


a compensating resistance (7), 
One pair of coils was 


few turns of wire, in series. 
wound from 33 gage manganin and the other from 
38 gage constantan wire. They were then impreg 
nated with glyptal cement and baked. To determin 
the Peltier coefficient, a pair of coils was connected 
in series with the Peltier junctions, so that the same 
current flowed through the coils and the junctions 
The coils were then inserted into the radio-balance 
one into each cup. The current was then adjusted 
in magnitude and direction to give approximate 
balance of the thermopile galvanometer, and the 
value of the current and residual galvanometer de- 
flection was determined. The coils were then inter 
changed between the two cups, the current reversed, 
and the new residual deflection measured. Knowing 
the value of the resistance of each coil, W’ and W”’ 
in eq (4) may be calculated, and d,, d,, and C have 
been measured directly. The value of s need only 
be known approximately as d,—d, is, in general, only 
a small correcting term. The value of P can thus 
be obtained directly with the help of eq (4). It 
should be noted that the energy dissipation of the 
coils appears in eq (4) as the difference between two 
wattages. The losses along the leads to the coils 
are, therefore, exactly compensated. 

The approximate value of s is determined by re- 
versing a suitably small current through the Peltier 
junctions connected in opposition, with no source or 
coil in either cup, and noting the change in gal- 
vanometer deflection. Ideally, dummy radium 
needles should be placed in the two cups. As WW" 
and W’”’ are now zero, s can be determined directly 
from eq (4), using the thermoelectric value of 7 
determined from eq (5). 

The values of P and r were determined in terms of 
the standard resistances, using the potentiomete: 
For this purpose, the lead wires were bared, and 
potential leads were attached at equal distances from 
the point of the coil, which had previously rested on 
the base of the radio-balance cup. As already pointed 
out, only W’—W’”’, and therefore R—r, is needed 
with accuracy. To obtain this, each coil was mounted 


on millimeter graph paper, and the distance from the 
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point of the coil to the point where the potential 
lends were soldered was determined to less than one- 
tenth of a millimeter. To enable this location to be 
made with accuracy and also to simulate the heat 
transfer of a radium needle, the platinum cores on 
which the two pairs of coils were wound, were turned 
to conical points, which protruded about a milli 
meter bevond the first turn of each coil at the end 
remote from the two current leads. 


4. Results 
Tarte l Resistances of calibrating coils 
Dista fror 
Cell " 1 end « Rp P 
) poter 
i Ohms Ohme Ohms 
{ ) 2 2H 0. 839 1. 421 
( tanta 2. 30 10.04 12 6. 822 
DD 7 9.103 2 24 6. 830 


The results obtained for ??—r for both the pair of 
manganin and the pair of constantan coils are shown 
in table | 

The results for the thermoelectric power (dkid7 
from which the Peltier coefficient (7dk/dT may be 
determined are shown in table 2. These can be com- 
pared with the value of 61.2 u~V/deg obtained at 
25° C by direct measurement of the thermoelectric 
power 


ranie 2 Vean thermoelectric power of Peltie inctions 


Residual | crhermoclectrs 
, ~— power, dE/d7 


Ohms " an pr deg 
1. 421 2 WY. 2 Ww 
1. 42) 25. 620 11.2 M4 

1. 421 26. 42 10.6 50.1 

1. 421 2s. i 45. ¢ 1 

6. 831 147 ‘ 09 

6. 831 ' ' 412 


The average of the first four values, for the man- 
ganin coils, is 59.12 107° wv/deg at 25.4° C; and the 
average of the last two values, for the constantan 
coils, is 59.10; 107° pv /deg at 24.3° C. An average 
value of 59.11 107° wy deg in the region of 25° C 
has therefore been assumed in caleulating the results 
that follow 

Results for the scale sensitivity, s, for different 
currents and for different cup fillings are given in 
table 3 

After these preliminary calibrations, little remained 
to be done but to use the radio-balance, and it trans- 
pired that it was not necessary to limit its use to only 
the conventional method of balancing the energy 
dissipation of a radioactive source against the 
Peltier cooling. Experiments have been carried out 
chiefly with four radium preparations in platinum- 
iridium needles. These preparations have en com- 
pared with the National Primary Radium Standard 
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and found, after correction for gamma-ray abserp- 
tion in the platinum-iridium walls, to have the fol- 
lowing values 


No. 3S015 1O.19 nw 
No. 359017 1.06 me 
No. BSO19 1.049 me 
No. 35920 O.USS mie 


If one of these sources be placed in cup A and a 
dummy needle in cup B, there will be an increase 
in the temperature of cup A and a corresponding 
deflection of the thermopile galvanometer. However, 
the largest of these sources, preparation 35915, 
causes a rise in temperature of the cup above its 
surroundings of only approximately 0.25 deg. C 
The loss of heat by radiation and gaseous conduction 
through the surrounding air will therefore be quite 
negligible, and the heat will be carried away mainly 
by conduction through the Peltier and thermopile 
leads. This quantity of heat will therefore be 
proportional to the temperature difference, or, 
conversely, the rise in temperature of the cup will 
be directly proportional to the quantity of energy 
dissipated within it. In other words, the galvanom- 
eter deflection should be directly proportional to 
the quantity of energy dissipated within it, and is 
therefore, in turn, directly proportional, in the case 
of radium, to the quantity of radium in the needle. 
It is necessary to qualify this statement in the case of 
preparations of radium that have been sealed at 
widely different dates and to take into account the 
small differences in energy dissipation caused by 
the growth of Po?”. With this qualification, however, 
the radio-balance has been found to provide an 
excellent method for quickly and accurately com- 
paring the ratios of quantities of radium, simply by 
observing the ratios of the galvanometer deflections, 

As an example of this method, the change in the 
galvanometer deflections on interchanging a radium 
source with a dummy between cup A and cup B is 
shown (fig. 5) plotted against milligrams of radium, 
for 3 of the 4 radium preparations used. Up to 5 mg 
of radium the scale deflection is a linear function of 
the mass of radium, but in order to preserve linearity 
up to 10 mg, it was necessary to convert scale 
deflections to angular deflections. 
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GALVANOMETER DEFLECTION. ON TRANSFER OF SOURCE, MINUTES OF ARC 


Figure 5 Galvanometer deflection in minutes of arc as a 
function of milligrams of radium transferred from cup A to 
cup B. 


Open circles are for sources 35019, 35017, and 35015. Solid circles are difference 
readings for source 35919 and source 35017 1.01 mg of Ra) and for source 35917 
and source 35015 («5.23 mg of Ra In order to reduce all readings to a com 
parable norm the milligram content has been increased by the ratio of column 4 
to column 5 in table 5 to allow for the fact that the radium content, as determined 
by gamma-ray measurements, does not take into account the growth of polonium 
The stated accuracy of the certification of the Radium Testing Laboratory is 
+0.7 percent his uncertainty in mass is approximately represented by the 
diameter of the open circles for the 10.1%mg source 10.32:+-0.07 meg, corrected for 
polonium 


It should perhaps be emphasized at this juncture 
that the energy emitted in the form of gamma 
radiation is largely lost from the radio-balance. 
With radium in equilibrium with its short-lived 
daughters, alpha radiation accounts for about 87 
percent of the energy emission, recoil nuclei for 
about 2 percent, primary beta radiation for about 4 
percent, and gamma radiation and secondary beta 
radiation for about 7 percent. Polonium in equi- 
librium would contribute a further 20 percent. 

The energies dissipated by sources 35915, 35917, 
35919, and 35920 were now measured by the transfer 
method between cups A and B, using the Peltier 
cooling and heating to equalize source and dummy, 
respectively. It later transpired that these sources, 
although purchased at the same time, had been 
sealed at different dates and it therefore became 
necessary to correct for the growth of the 5.3-Mev 
alpha particle due to Po”. This was done by using 
a relation derived by Curie and Yovanovitch [4], 
who neglected the effect of the short-lived intermedi- 
aries between radium and radium D and between 
radium D and polonium. This relation, the numeri- 
cal constants of which were later revised by Saniele- 
vici [5], Zlotowski [6], and then by Jordan [7], using 
more accurate half-life data, gives Q.,,, the heat 
generated by the Po*" alpha particles, as a function 
of the time ¢ in years as follows: 


1.0175¢e7-%"!+.0.0173e""™"" eal go 
hr=' (6) 


Qa ,=27.26(1 


| 
| 
| 








The results obtained are summarized in table 
All sources are corrected for polonium growth {) 
October 1953. 

The ratios of these sources obtained by gamma-r 
comparisons in the Bureau's Radium Testing Labo. 


and by the radio-balance are summarized 
The stated accuracy of the Radium Testi: 
t0.7 percent, 

Values obtained at other 
emission of sources 35917 and 
table 6. 

Values for source 35917 were also determined 
using cup A and cup B individually. For this pur- 
pose, it is necessary to determine the resistance 7, 
and Ry, of the Peltier junctions to cups A and B 
Using eq (2), these were found to be 46.30 10~ and 
46.07 <10~° ohm, respectively. As has been men- 
tioned previously, however, it is difficult to remove 
completely the residual thermopile zero deflections 
for each cup individually. In this determination 
therefore these zero deflections were determined 
initially, and the values for the source energy emis 
sion were reduced to these afterward by extrapolation. 
However, this method is tedious compared with the 
transfer method. The results obtained by using eq 
(1) are shown in table 7. 

One final point of note is that the radio-balance 
can be used to compare two nearly equal sources by 
the transfer method without using any Peltier current 
atall. From figure 5 the scale sensitivity for transfer 
of a mass of radium from cup A to cup B is 0.0193 
mg/mm. On placing source 35919 in cup A and 
source 35920 in cup B, an out-of-balance deflection of 

1.9 mm was observed. On interchanging the 
sources, an out-of-balance deflection of —1.4 mm 
was observed, giving a total change on transfer of 3.3 
mm. In actual practice, twice this deflection is 
observed on reversing the current through the gal- 
vanometer, but for convenience in recognizing 
whether cup A is warmer or cooler than cup B, these 
are reduced to half the deflection to the right (++) or 
left (—) of zero. 

A change of deflection on transfer corresponds 
therefore to a difference between sources 35919 and 
35920 of 0.064 mg, with 35919 the greater (a positive 
deflection corresponds to cup A, being at a higher: 
temperature than cup B). 

Thus if we assume that source 35919 consists of 
1.049 mg of radium, as certified by the Radium Test- 
ing Laboratory, then source 35920 comprises 0.985 
mg (compared with 0.988 mg as certified by the 
Radium Testing Laboratory). 

This application is also illustrated for widely 
differing sources in figure 5, where points derived by 
taking differences between sources 35915, 35917, and 
35919 have been plotted. 

A series of experiments was also carried out with 
the radio-balance, using individual cups in an evac- 
uated bell jar at approximately 5x10~-° mm Hg 
pressure, with a view to determining whether the loss 
of heat from the cups by radiation and gaseous con- 
duction was of significance. In general, it was found 
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R515 April 6, 1951 10. 19 42.674 118 141.1 151. 2 140.3 
5017 November 6, 1952 1 21. S87 1.1 Ta 0 151.0 15 
er 35919 December 22, 1049 1.049 + 50 a0 159.9 182.4 140.4 
© sewn December 22, 1040 0. ORS +. SAS 52 151.7 153 1™) 
‘nm in 
* Average energy dissipation per milligram of radium (less P,.™° contribution) equals 149.9 gwime +0.4 percent, which equal to 128.9 cal ¢ w rt 
d by percent + 0.4 percent represents a 95-perent confidence limit 
pur- = : , oud 
oR Paste 5. Measured relative values of radium sources pump oil originated, however, from the experiments 
, 
dB with the apparatus In a vacuum, but its use was 
and NBS Radium Testing Ce a er ne continued because of the obvious advantage of a low- 
~worator . . . . . 
men icioantias vapor-pressure liquid from the point of view of avoid- 
nove Soures ices ing decreases in temperature due to evaporation. 
tions Value Ratios Pa ao Ratios ° > 
ation contribution 5. Comparison With the Mound Laboratory 
ned me - és Calorimeters 
Mts S591 10.19 1541.1 1521.8 1 
$5017 406 0. 487 749.0 746.0 0. 490 ° . 
tion s5019 1. 049 1020 159.9 156.8 1030 As the total energy emission of two radium needles, 
1 the 5920 0. 968 eure ASI. 3 148.6 6 | | the milligram content of which had been measured 
ig eq in terms of the National Primary Radium Standard, 
Taste 6. Energy dissipation of two radium sources had recently been measured in the calorimeters of 
lance the Mound Laboratory, it has been possible to com- 
s by , pare their results with those obtained with the radio- 
energy dissi 
rrent me atten tame Date of balance. 
nsfet ee The total radium content of these two radium 
0193 sources as measured by the Radium Testing Labora- 
and a tory was 10.20 mg. 
on of eer ag Spent, ae The Mound Laboratory obtained results both 
the 45017 748.7 September 24, 1953 | with their microcalorimeter and their macrocalo- 
O17 ‘7 ~ptember 25, 105 . Pe 
mm 1 46 Once Be | rimeter. With the former they found the total 
of 3.3 aie a sasha. Shi energy output of the two needles to be 1.634 uw, 
wit 159.0 Septem he , 1s > . ° . . 
ym is 35019 158. § Do with 0.1 em equivalent lead absorption; with the 
501 ¢ LS. 6 Do ; “4 ; 
_gal- oaeae oan 8 October 22, 1988 latter, they found it to be 1.65 ww, with 1.1 em 
izing equivalent of lead absorption. 
these * The mean value for source 35917 is 746.3 uw, with a Initially, it was assumed by the Mound Labora- 
\ standard de ion of ( vercent. he mes © fo . . 
) or many Pp =n ye Sp Se ae | tory that these needles had been sealed recently, 
0.8 percent However, it has transpired that they had been first 
* No dummy used in other cup r . 7% 
onds tested by the National Bureau of Standards in June 
and Prarie 7 Determination of the enerqu emission froma radium 1939, and that considerable correction had to be 
Rg source, using the radio-balance cups individually a a >, 210 
itive made for the growth of Po*®. 
het By using a value of 4.185 j/cal and the value of 
> . ; . 
Balancing Residual Energy 10.20 mg of radium, the Mound Laboratory results 
Cup eurrent galvanometer Temperature dissipation,» daa et 137 Se | per | —w "ee om | rit ‘ter. 0.1 
ts of leflection * PC-—CU |} reduces oO Od. cal £ iT mhicrocatorimeter, VU. 
est- / em Pb equivalent) and 139.2 cal g™' hr-'(macro- 
) ORS . J. “er * a calorimeter, 1.1 em Pb equivalent), both uncorrected 
Oe 7. 76 +4) 2 735. 2 | > . ° 
i | \ Ws 325 11 207.8 742.1 for Po*®® growth and therefore including all energy 
Af 0 007.7 re rs : 4 
. oa » or gd contributions down to that due to the Po*"® alpha 
del particles. 
aery ° ° 
' : * The initial “zero” deflections for cups A and B, with a dummy 5-ing The cups of the radio-balance are 0.04 em thick: 
d by needle that had been left to come into thermal equilibrium in each, wer: . 7 _ . -“ . 
ond 1.3 and —1.0 mm, respectively. Thus the value of source 35917 as deter. | in each case the wall thickness of the platinum- 
mined in cup A is, by extrapolation, 743.2 ¢w. The value for cup B is 752.0 | Senals _ . . Rw > . —— “or ; 
ww. The average value for cups A and B is 747.6 aw (or 745.4 uw after correct- iridium nee dle was 0.05 em. a As the mass absor ption 
J ing for the growth of polonium) compared with 746.5 #w obtained by the coefficients of copper, platinum, and lead are, for 
with transfer method (also corrected for polonium growth).f : w vin heornti 
the purpose of making a gamma-ray absorption 
vile that the radio-balance was much more unruly, the | correction in this energy range, sufficiently near to 
g extraneous thermal effects being more difficult to | equality, and as the needles consist of 90 percent of 
}s085 track down and eliminate. This method of using the | platinum, it will be sufficiently accurate to reduce 
ee radio-balance was therefore abandoned in favor of | the absorption of the platinum-iridium needles and 
yun | the methods described. The use of silicone vacuum- | the copper cups for gamma radiation to their equiva- 
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lent total absorption in lead using the densities of 
copper, platinum, and lead. This gives a gamma- 
ray absorption of the needles and cups equivalent 
to that of 0.126 em of lead 
Assuming that the radium needles had been sealed 
for 14 vears prior to their measurement by the 
Mound Laboratory between December 1952) and 
February 1953, the contripution due to Po*® would 
have been 9.5 cal g°' hr-', as caleulated from eq (6) 
We have, therefore, the following results for the 
energy output of 1 g of radium and its products 
down to radium D 
Mound microcalorimeter with 0.1 
em Pb equivalent 
Mound macrocalorimeter with 1.1 
cm Pb equivalent 
128.0 cal go! he NBS radio-balance with 0.126 em 
Pb equivalent 


128.3 cal go! be 


120.7 cal go! he 


These results have been reduced, for the purpose 
of comparison, to the same value of equivalent 
absorption, using the measurements of Zlotowski |6), 
who determined the energy emission of four radium 
sources calorimetrically by using different thicknesses 
of lead absorber around the sources. Zlotowski's 
results are shown in figure 6, where the heat absorbed 
in the calorimeter is plotted as a function of thickness 
of lead absorber. By using this curve the Mound 
Laboratory results may be reduced to an equivalent 
thickness of lead equal to 0.126 em and compared 
directly with the NBS value. The results then 
obtained are 
128.5 cal wt he 
125.4 cal g@' hbr—'! 

128.9 cal g-' hr—-'! 


Mound microcalorimeter 
Mound macrocalorimeter 
NBS radio-balance 

The probable error of the mean value obtained 
from 10 readings with each of the Mound calorim- 
eters is +0.5 percent. The probable error of the 
mean value obtained from 4 readings of the NBS 


radio-balance is + 0.2 percent 














140 
© = 
+8 
P 
wal / 
- “ 
o 
Vv 
o 135 O ; ; 
@ ; 
a 7 
° P 
” 
@ é | 
< / 
f | 

a p 
we ; 
= ° 

’ 

130 | 
° 2 4 6 a 
EQUIVALENT LEAD ABSORPTION, cm 
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function of equivalent absorption in lead, after Zlotowski [6] 


The point at 7 em of 199.6 cal ¢' br-' is the extrapolated value for the total 
heat emission of radium-226 and its daughters down to Ra C 


6. Conclusion 


It is proposed to apply the radio-balance described 
in this paper to the measurement of the average 





energy of beta emitters and also to a number of oth: 


appropriate problems. A new radio-balance is bein 
constructed with a view to comparing the U.S. an 
Canadian National radium standards. Two sue} 
standards were compared with quite  promisin 
results in the radio-balance here described, but th 
cups of this radio-balance are somewhat too sma 
for these sources. It is interesting to note, howeve: 
that the radio-balance is so responsive that a con 
parison agreeing well with the known data was ca) 
ried out between the U.S. and Canadian Nationa 
radium standards in as short a time as | hr 

From table 3 the average value of the seale sen 
sitivity for transfer is 0.332 mm uw. As deflection 
may be read to about 0.1 mm (i. e&., 0.2 mm or; 
reversal of the current through the galvanometer 
it should be possible to determine an energy emission 
to about a third of a microwatt. Thus a 30-uw 
source should be susceptible of measurement with 
an accuracy of about | or 2 percent. Such a valu 
corresponding to about 0.2 mg of radium, probably 
represents the lower limit of usefulness of the radio 
balance 
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Some Factors Affecting the Dimensional Stability of the 
Silver-Tin-(Copper-Zinc) Amalgams* 


J. A. Mitchell,' I. C. Schoonover, George Dickson, and H. C. Vacher 


\n investigation of the X-ray diffraction characteristies and of the dimensional changes 
of silver-tin-(copper-zine) amalgams during and after their hardening period indicated that 
uncombined mercury Was present in the amalgams after the initial solidification and that 
this uncombined mercury disappeared during the period when greatest dimensional changes 
occur On the basis of the data obtained, it is postulated that the expansion during the 
hardening of amalgam results from diffusion of uncombined mereury throughout the material 
the subsequent shrinkage then results from combination of this mercury with existing phases 


or with residual alloy. 


1. Introduction 


Although an enormous amount of experimental 
work has been done on the silver-tin-(copper-zine) 
amalgams [1,2,3,4,5,6] no mechanism has been 
proposed that will satisfactorily explain the initial 
contraction, subsequent expansion, and final con- 
traction that occur during the first 24 br after 
amalgamation. Various investigators have = at- 
tempted to explain these changes in terms of mercury 
content of the amalgam system. Their difficulty 
may have resulted from the failure to differentiate 
between total mercury and combined mercury in this 
system. It is probable that in the hardening of the 
heterogeneous mass not all the mercury is in the com- 
bined form. Chemical analysis of the amalgam 
would not distinguish between combined and = un- 
combined  mereury. Furthermore, —photomicro- 
graphic evidence [7] indicates that the rate of attack 
on the alloy particles by mercury may vary widely, 
depending upon the conditions under which the speci- 
men was prepared. The possible diffusion of this 
mereury, if it is present, and its subsequent reactions 
with the alloy may account for the dimensional 
changes observed in the hardened amalgam. 


With this in mind, the investigation was under- 
taken to determine (1) whether uncombined mercury 
was present after the amalgam had hardened, (2) 
whether the uncombined mercury, if present, con- 
tinued to take part in any reaction, and (3) whether 
the dimensional changes that oecur in hardened 
amalgam could be correlated with this reaction. 


2. Experimental Approach and Procedures 


X-ray diffraction methods cannot be used to detect 
the presence of free or uncombined mercury in an 
amalgam at room temperature because the mercury, 


“A dissertation submritted to the faculty ef the Graduate School, Georgetowr 
University, in partial fulfillment of the requirements for the degree of Master 
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(iuest worker at the National Bureau of Standards, and member of Dental 
orps, Department of the Navy 
he chemical symbols Ag-Sn-(Cu-Zn) are used hereafter. Copper-zine is 
placed in parentheses because these elements are not always present in the above 
thentioned amalgams 
Figures in brackets indicate the literature references at the end of this paper 


being a liquid at this temperature, does not have a 
crystalline structure and does not give an X-ray 
diffraction pattern. However, if the mereury ts 
chilled below its solidification point, typical lines 
representing mercury are detected by diffraction 
methods. Therefore, if uncombined mereury is 
present in amalgam immediately after hardening, it 
should be possible to determine its presence by X-ray 
investigation of amalgam specimens held below the 
freezing temperature of mercury. Furthermore, it 
should be possible to study the rate of disappearance 
of the uncombined mercury with time by periodic 
examination of the specimen at temperatures below 
the freezing point of mercury. 

X-ray diffraction patterns of the amalgams at low 
temperatures were obtained from specimens that 
were manually condensed, using standard dental 
procedures, in a stainless-steel mold and specimen- 
holder combination (fig. 1). Immediately after con- 
densation, the amalgam specimens were frozen by 
introducing liquid nitrogen into the chamber beneath 
the specimen. An iron-constantan thermocouple 
was inserted in the specimen holder, near the speci- 
men for temperature determination 

Three commercial amalgam alloys were used in the 
investigation. Amalgams were produced by tritu- 
rating the alloys with mercury in a mechanical 
amalgamator, using the mercury-alloy ratio of 2 to 3 
for alloys | and 2, and a ratio of 2 to 3.6 for alloy 3. 

A high-angle spectrometer goniometer [8] was used 
to obtain chart records of the copper K-alpha diffrac- 
tion. The diffraction unit was operated at 35 kv and 
a tube current of 20 ma. The angular speed of the 
scanning arm was | deg/min, and the chart speed was 
‘ mn./min, 

X-ray diffraction characteristics were obtained 
with the specimen and its holder maintained at the 
constant reduced temperature (—125° C) during the 
making of chart records in a 26 range of 80 to 28 deg 
The specimen was then warmed to mouth tempera- 
ture (37° C) and aged at that temperature until the 
next low-temperature diffraction-chart recording was 
made. 

The X-ray diffraction technics described above 


| were used also to determine the effect of addition of 
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SECTION A,, A, 
Fieure |! tmalgam mold and specimen holder for low- 
tem perature X -ray diffraction studies 
Holder enclosed in a cyl ndrical box. Specimen retained in holder and alter 
nately cooled to —125° C by liquid nitrogen and warmed to 37° C by a wire 
wound resister regulated with a variable autotransformer 
A, Liquid nitrogen intake; B, liquid-nitrogen chamber; C, jets; D, specimer 
FE, thermocouple rest; F, arm fastened to goniometer shaft 
a 
8 
r 
Fiavrr 2 Measurement of dimensional change of fully 


hardened amalgam specimens after addition of mercury 


Mercury added internally and externally 


A, Dental interferometer; B, specimen; C, mercury; D, iron cup. 


mercury to fully hardened amalgam specimens. 
Seven-day-old specimens were immersed in mereury 
for 12 hr, removed and cooled to 125° C, and 
examined by X-ray diffraction. The specimens were 
kept at 37 °'C during the next 24 hr, exe ept for peri- 
odic cooling to — 125° C while the diffraction charts 
were made. 

To examine the effect of uncombined mercury on 
dimensional changes of hardened amalgam, mercury 
was added both externally and internally to hardened 
amalgam specimens, and their dimensional changes 
were measured by means of dental interferometers 
[9], (fig. 2). 








iron cups were positioned in interferometers. Me: 
cury was added to the cups, and the dimensions 
changes were observed over a period of 15 days o 
two series of specimens held at 21° and 37° (¢ 
respectively. In other experiments a small hole wa 
drilled in the specimen and mercury injected int 
the hole in the amalgam by means of a hypodermi 
syringe after the specimens had been positioned 

the interferometers. Dimensional changes were ther 
observed over a period of 15 days. Temperatur 
were maintained by constant-temperature air baths 

In two of the above experiments the supply of 
mercury was removed from contact with the speci 
mens after the period of 15 days, and dimensional 
changes were observed for an additional period of 
26 days. 

To obtain information on the effect of temperature 
on the reactions occurring during the hardening of 
the dental amalgams, dimensional changes of amal- 
gams prepared and observed at various temperatures 
were also determined interferometrically. 


3. Results 


3.1. X-ray Diffraction Examination of the 
Ag-Sn-(Cu-Zn) Amalgams 


As shown in table 1, one of the three commercial 
amalgam alloys used in the investigation was a 
silver-tin-zine alloy, one a silver-tin-copper-zine, and 
the other a silver-tin-copper alloy. Upon examining 
their X-ray diffraction patterns (fig. 3), it is evident 
that the variation in copper and zine content has 
little effect on the diffraction patterns as the position 
of the principal peaks are identical in all three charts 

TABLE I. 


Chemical composition of alloys investigated 


Alloys 
Elements 
l 2 $ 
Ag 74 70.2 70.1 
Sn 25.0 24.6 26.5 
Cu 0.0 24 4 
Zn o.5 18 0.0 


X-ray diffraction charts of pure mercury at — 125 
C show a number of well-defined peaks (fig. 4, a). 
A specimen of amalgam containing the Ag-Sn-Cu- 
Zn alloy was made at room temperature (25° C 
nese: in position, and its X-ray pattern determined 
at 25° C 5 min after the start of mix (fig 4 (b 
The ¢ as ‘teristic peaks of the uncombined mercury 
can be seen as represented by the solid-black areas. 
As the specimen is aged through longer periods of 
time at 37° C, the characteristic peaks of the uncom- 
bined mercury tend to disappear, and after 6 hr they 
are scarcely detectable (fig. 4 (c), and fig. 5). This 
time period corresponds generally to the period 
during which significant expansion and contraction 
are commonly observed. Simultaneously with this 


Hardened specimens of amalgam in | reduction in uncombined mere ury, there is a disap- 


186 

















ral 
ires 


cial 
Ss ii 
and 
ing 
lent 
has 
Lion 
rts 


C), 
ned 
(b 

ury 
eas. 
s of 
om- 
hey 
This 
riod 
tion 
this 
sap- 











\ 








80 


60 


40 


20 


80 


60 


40 


20 


60 


40 


20 


The similarity of the charts for the three alloys indicates that none of the lines shown are dependent upon the copper o 
alloys 
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N-ray diffraction charts (with a scale factor of 4 


Figure 3. made at 25° ( 


r zinc content of the 
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Ag-Sn-Cu-Zn AMALGAM 
AGED 5 MINUTES AT 25°C 
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Fiavuar 4 \-ray diffraction charts made at 126° C. 


ur\ b and ¢, Scale factor of 4 used for the remainder of the procedures. Height of the peaks o 


Seale factor of 32 used for obtaining the tracings of m 
y diffraction charts made at 37°C. Lines representing ice formation are labeled “ice 
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X-ray diffraction charts made at — 126° C, 


Fieure 5 


(See legend for fig. 4, b and « 
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Ficure 6, X-ray diffraction charts made at — 125° C. 


(See legend for fig. 4, b and c) 

















pearance of the lines representing the original alloy 
y)' and an increase in the 7; and y2 phases (fig. 4 (b) 
and (ce), and fig. 5). This indicates that the disap- 
pearance of the free mercury lines was through 
combination of the mercury with the alloy to form 
mercury compounds. 

X-ray diffraction charts made periodically at re- 
duced temperatures after the addition of mercury to 
fully hardened amalgam showed intense mercury 
lines in the specimens observed immediately after 
removal from the mercury bath (fig. 6). These 
mercury lines were not detected after the specimens 
had been held at 37° C for 24 hr or at 21° C for 3 
davs. In figure 6 it can be seen that mercury on the 
surface of the specimens reduces the intensity of the 
lines representing the original alloy and y, and y, 
phases. After aging for 24 hr, the mercury has com- 
bined with the original alloy as indicated by the 
absence of y and uncombined mercury lines. 

Consideration was given to the possibility that the 
intermediate phases in the Ag-Sn-(Cu-Zn) amalgams 
might have a higher mercury solubility at room tem- 
perature than at temperatures below the freezing 
point of mercury. Under such conditions, crystalline 
mercury at low temperatures might be formed by 
rejection from a solid solution existing at room tem- 
perature. However, the disappearance of the lines 
representing uncombined mercury indicates that such 
a rejection does not occur to a significant extent if 
at all. 


3.2. Dimensional Change Determinations of the 
Ag-Sn-(Cu-Zn) Amalgams 


The addition of mercury to the hardened specimens 
produced large expansions, whether the mercury was 
added internally or externally (fig. 7). Specimens 
maintained at 37° C expanded much more than those 
maintained at 21°C. The great difference in expan- 
sion between the two specimens placed in mercury 
baths at 37° C may result from imperfections or other 
variations in the individual specimens that affect the 


rate of diffusion of mercury into the amalgam. Only 
very slight dimensional changes occurred after 
removal of the supply of excess mercury from 


hardened specimens (fig. 8). 

Amalgam specimens prepared and maintained at 
either 21° or 37° C showed a rapid initial contraction, 
followed by some expansion during the first 24-hr 
period. Average values for dimensional changes are 
plotted against time in figure 9, 10, and 11. 
Ag-Sn-Zn amalgams prepared at 21° C and main- 
tained at 37° C gave an over-all shrinkage, whereas 
those prepared and maintained at 21° C gave an 
over-all expansion (fig. 9 and table 2). Ag-Sn-Cu-Zn 
amalgams prepared and maintained at 37° C also had 
greater initial contraction, lower maximum expansion 
and greater over-all contraction than those prepared 
and maintained at 21 °C (fig. 10 and table 2). 
Ag-Sn-Cu-Zn amalgams initially observed at 21° C 
and then subjected to a temperature increase to 
C shrank more than those maintained at 21° C 
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(alloy 2) resulting frem addition and removal of mercury 


0 days, mercury added internally through hole drilled into specimen; 15 days, 
excess mercury removed from hole drilled into specimen 
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Ficure 9. Comparison of dimensional changes of Ag-Sn-7n 
amalgam (alloy 1) at 21° C and at 37° C 
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erence 


during the observation period (fig. 11 and table 2), 
Specimens held at 37° C and then subjected to a 
temperature decrease to 21° C shrank less than those 
maintained at 37° C (fig. 11 and table 2). Ag-Sn-Cu 
amalgams observed at 37° C also shrank more than 
those observed at 21° C (table 2). 

Averages of the differences in dimensional changes 
of the Ag-Sn(Cu-Zn) amalgams at 21° and 37° C are 
given in table 2, together with their standard errors. 

The standard error in each case depends on the 
number of specimens and on the precision of a single 
measurement. The latter was evaluated separately 
for each period of hardening, by pooling the variances, 
obtained for all amalgams, of the measurements 
made at that setting time 

In table 2 it is also shown that the fiducial ob- 
servation (zero reading) taken 10 min after start 
of mix yielded a significantly lower expansion or 
greater contraction value over a 24-hr period than 
the 15-min fiducial observation. This variation 
results from the fact that a larger portion of the 
initial shrinkage is observed when the 10-min 
fiducial reading is used 


4. Discussion of Results 


The observations of dimensional changes, along 
with the X-ray diffraction evidence that uncombined 
mercury is present in the amalgam, tend to sub- 
stantiate the theory that the mechanism of dimen- 
sional change is associated with the presence of the 
uncombined mercury 

The following mechanism is postulated to ex- 
plain the physical phenomena observed. The initial 
shrinkage always observed results from the formation 
of compounds of mercury and alloy having a smaller 
volume than the sum of the original volumes of the 
alloy and mercury. The subsequent expansion that 
is observed results from the diffusion into new areas 
of uncombined mercury entrapped in spaces through- 
out the amalgam. The reaction of this uncombined 
mercury with existing phases or with residual alloy 
to form compounds of lesser volume would then 
account for the contraction that follows the expansion 
of the amalgam. 

This mechanism would account for the dimensional 
changes that occurred on the addition of mercury to 
hardened amalgam. When mercury is added to 
hardened amalgam the diffusion of this mercury 
throughout the specimen would produce the large 
expansion observed. The slight shrinkage ob- 
served immediately after the removal of the supply 
of mereury could result from the combination of 
free mercury in the specimen with existing phases or 
with residual alloy. 


The greater shrinkage of amalgam mixed and 
maintained at 37° C rather than at 21° C may result 
from a more rapid chemical combination and dis- 
appearance of the uncombined mercury at elevated 
temperatures 

The proposed mechanism would also account for 
the shrinkage of amalgam that occurs from over- 
trituration or excessive working during condensa- 
tion, in that either of these treatments would tend 
to reduce the amount of uncombined mercury by 
either furnishing new surfaces for reaction with the 
mereury or by actual mechanical removal of the 
mereury These treatments would also tend to 
distribute the uncombined mercury present more 
uniformly throughout the amalgam so that there 
would be fewer mercury-rich areas from which dil- 
usion could take place. On the basis of this ex- 
planation, a reduction in size of the alloy particle 
or a heat treatment, which would make the particles 
more reactive, would tend to cause a reduced ex 
pansion or a shrinkage of the amalgam 


5. Summary 


Residual, or uncombined, mercury has been 
demonstrated by X-ray diffraction studies to be 
present in amalgams after their initial solidification 

This uncombined mercury disappears during the 
period when greatest dimensional changes occur in 
the hardening of the amalgam system 

A possible mechanism for the dimensional changes, 
based on the presence and disappearance of the 
uncombined mercury, is presented 
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A Colorimeter for Pyrotechnic Smokes* 


Isadore Nimeroff and Samuel W. Wilson 


\ tristimulus photoelectric colorimeter has been designed to measure the chromaticities 


of pyrotechnic smokes having highly saturated colors 

designed as a 
CIE 
instrument. 
little difficulty was encountered in designing filters to approximate 


wall reflection, the smoke chamber was 
designed to approximate the daylight 
filter-phototube combinations of the 
modal, 


To minimize edge effects and rear- 
36-inch cube. Three filters were 


standard observer functions for the source- 
As the CIE 7- and 2-functions are uni- 
these functions. 


The CIE a-function, which is bimodal, required a divided filter to approximate its two lobes. 


The filter constructed was composed of two sectors, one 
7-function. 
, using Munsell papers as standards 


other approximating the “red lobe’’ of the 


tained with this colorimeter 
1. Introduction 


In 1931 the International Commission on Illumi- 
nation [1]' (CLE) adopted a system for reducing 
spectrophotometric data to color coordinates in 
three-space. The reduction of such data for a single 
specimen entails a considerable amount of computa- 
tion. To avoid this computation, Twyman and 
Perry [2] suggested in 1930 that three filters dupli- 
cating the then proposed standard-observer func- 
tions could be used in a photoelectric instrument. 
Several such instruments, now called photoelectric 
colorimeters [3], have been designed since then. 
The photosensitive elements in colorimeters have 
been barrier photocells [4, 5], vacuum phototubes [6], 
and multiplier phototubes [7, 8] 

This paper describes a colorimeter, using multi- 
plier phototubes, deslanad at the request of repre- 
sentatives of the Army Chemical Center to be used 
for evaluating the colors of pyrotechnic smokes. 
This instrument was needed to determine progress 
in the development of saturated smoke colors and to 
check on reproduction of specified smoke colors. 


2. Filter Design 


The colorimetric spectral functions approximated 
in this colorimeter by combined source-filter-photo- 
tube response are those established by the CIE. 
The theory of the design of filters is based on two 
propositions (1) that the three CIE observer func- 
tions represent the colorimetric system of the average 
normal observer, and (2) that each of these functions 
can be approximated by the proper selection of 
optical filters in combination with a light source and 
a phototube. As the instrument is intended to 
evaluate the color of nonfluorescent smokes as if 
they were viewed in daylight (6,500° K), even though 
in the colorimeter they are illuminated by incandes- 
cent-lamp light (2,780° K), the proper correction was 
incorporated in the filter design: Had the colorim- 
eter been intended to evaluate the color of fluores- 


—— 


*This project was done under contract with the Chemical Corps, Department 
of the Army. 
' Figures in brackets indicate the literature references at the end of this paper. 





*, the 


been ob- 


approximating the ‘“‘blue lobe’ 


Satisfactory results have 


cent, as well as nonfluorescent, smokes as viewed 
under daylight, the approximation of the CIE func- 
tion would be accomplished by phototube and filter 
combination alone, but the illuminator would then 
be required to yield artificial daylight. 
The CIE spectral-response functions, 7, 


y, and = 


under daylight spectral conditions, /,, are respec- 
tively, 7E,, YE,, and ZE,. If the illuminator yields 


spectral distribution, /,, the spectral responses, S,, 
S,, and S,, of the phototubes for each of these func- 
tions must be altered by filters having spectral trans- 
mittances 7,, 7,, and T,, respectively, in order to 
refer to viewing by daylight. This condition may 
be stated mathematically: 


E,T.S,K, 
E,T,S,K, 
E,T.S.K, 


VE, (1) 
2E, ! 

where K,, A,, and K, are proportionality constants, 
independent of wavelength for their respective func- 
tions. The desired relative spectral transmittances 
then are determined by solving eq (1) for K,T,, 
K,T,, and K,T,. 

The filter design is accomplished by determining 
the proper combination of colored glasses. Where 
the transmittance, 7), of a glass of stock thickness 
dy, is to be changed to a transmittance, 7, the fol- 
lowing equation must be solved for the new thick- 
ness, d;: 

T= (To) k?)'- 4 


0, (2) 
where *& is unity minus the Fresnel reflectance of the 
glass. The approximate value of & for glass is 0.96. 
To minimize the work entailed in this computation 
a thickness-transmittance nomograph was used [9]. 

Figure 1 shows the spectral-response curves of 
several multiplier phototubes (RCA 5819) as cali- 
brated by the NBS Radiometry Laboratory. The 
first phototube purchased had the response ‘shown 
by curve A. Because of the relatively high response 
of this phototube in the region about 580 mu the 
spectral transmittances of required filters for CIE 
functions must have absorption bands in this region. 
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phototubes 


To duplicate these complex transmittance curves 
the technique of divided-disk filters was devised. 
Fortunately, the phototubes purchased later had 
normal responses, as shown by curves 3, 5, and 12, 
of figure 1, thereby simplifying the filter-design 
problem. The phototubes having these curves are 
the ones used in the colorimeter for functions Y, 2, 
and 2, respectively. 

As no combination of glasses has ever achieved 
duplication of the bimodal Z-function in a single 
filter, the divided-disk technique devised for the 
complex transmittances required for phototube A 
was applied to the design of the bimodal 7-function. 
This division requires that one part of the single filter 
approximate the “blue lobe”’ and the other part 
approximate the “red lobe” of the Z-function. The 
transmittance of the ideal filter, 7, for the Z7-function 
may be divided into its two lobes, 7,, and T7,,, thus 
T.=T,4+-T.,. In terms of relative transmittance, 
this becomes 

K,T, (3) 


KT 5+ T1). 


The actual approximation consists of a disk, | 
area A, with one fraction, a/A, of the total area 
having spectral transmittance 7.,, and the remaining 


fraction 1—(a/A), having spectral transmittance | 
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wes Then the transmittance, 7/7, of the act 
filter is 
wa Dae ¢ eee ow 
1 T".4 _ A { 
- A : 
The conditions that the spectral lobes of th 


approximate filter be properly proportional to thos 
of the ideal filter are: The summation under eac! 
lobe of the approximation must equal the correspond 
ing summation under the ideal filter, thus 


TS K.T.=DK.1 
ST K.T.=SK.7 


The ratio of the summation under the two lobes i 
the approximate filter must equal the ratio of th 
summation under the two lobes of the ideal filter 
thus 


~ @ ww 

a4! S7 
; fh) 

wv A a ‘a S357 


But substituting the values of 7), and 7 from eq 


(5) into eq (6) gives 


Ky | . 
A kK,+Kn 


The area of the disk of radius r may be divided 
into either segments or sectors. Consider the seg- 
ment division. The total area, A, is zr’, and the 
area of the segment, a, is (r*/2)(@—sin @), where @ 
in radians, is the central angle subtended by the seg 


ment, as shown in figure 2,a. Then 
a ‘ 
(@—sin @). &) 
i 2 


After the value of @ is computed, the segment is 
determined by computing, r=r cos(@/2). 

Although relatively easy to construct, a filter with 
this type of division of the disk has two major diffi- 
culties. The disk, once constructed, is useful only 
with a circular opening of radius r. As the segment 
is displaced from the center of the disk, one must 
assume that the spectral response of the photosensi- 
tive surface of the phototube is completely uniform. 

Consider now the sector division of a disk filter 
The area of the sector (see fig. 2,b) is equal to r’@/2, 
where @, the sector angle, is in radians. Then 
a/A=6/2. If the sector angle, expressed in degrees 
is p, then 
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A 360° 
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a b 
Figure 2. Divided disks for t-function filt 


a, Segment division; b, sector division 


This type of division, although more difficult: to 
construct, does overcome the major difficulties of 
the segment division. The spectral character of the 
transmitted light is independent of the radius of the 
disk. One must assume only that the spectral 
response of the photosensitive surface of the photo- 
tube is radially symmetrical. This is a reasonable 
assumption for the head-on-type phototube used, as 
the photosensitive surface is evaporated onto the 
center of the cathode, thereby vielding a close ap- 
proach to radially symmetrical sensitivity. This 
radial symmetry can readily be tested, furthermore, 
by rotating the sectored filter in its own plane and 
observing if the output of the phototube changes 
For either type of division the center of the disk 
must be carefully positioned over the center of the 
photosensitive surface 

Because of the advantages outlined above, the 
sector division of the filter was used. On the basis 
of the application of eq (2) to the published spectral 
transmittance of stock glass, the combinations vield- 
ing the best expected approximation to the two lobes 
of the 7-function were chosen. The design data are 





RELATIVE TRANSMIT TANCE 





\. 
£ i Wy Wa 
400 450 $00 550 600 650 
d» (my) 





Jj \ 4 f , + 
700 


Fieurre 3 Spectral transmittance of the desired (@), expected 
A), and actual filters for CITE F-function 


shown in table 1. The sector angle was determined 
by substituting eq (7) into eq (9), giving 


p — 360 K,,/(K.,+Ke) (10) 


for the Z-function filter designed for this colorimeter, 
K,,=0.915, K,,=5.483, and p=51.5°. The expected 
relative transmittance, A777, for this filter is shown 
in figure 3 compared with the desired relative trans- 
mittance, A,7,, and the achieved relative trans- 
mittance, A,7?, of the actual filter. The actual filter 
proved to be a better approximation to the desired 
filter than had been expected from the design data, 


TABLE 1 r-Junction filter design data 
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As the CIE 7%- and 7-functions are unimodel 
little difficulty was encountered in designing filters 
to approximate these functions. For the source- 
filter-phototube combinations of this colorimeter, the 
approximations (light lines) to the CIE tristimulus 
functions (heavy lines) are shown in figure 4 


3. Description of Apparatus 
3.1. Chamber 


Most of the colorimeters heretofore designed have 
been used for measurement of colors of opaque or 
transparent media. For these media, where edge 
effect is negligible, the size of the illuminated and 
viewed area of the medium was determined by the 
sensitivity of the instrument. The instrument here 
described is being used to measure the color of trans- 
lucent media, pyrotechnic smokes. As the smokes 
are translucent, that is, light scattering as well as 
transmitting, and reflecting, a rather large chamber 
having a uniformly illuminated window is required 
to minimize edge effects. 

Figure 5 shows a general view of the apparatus. 
The chamber is a l-vd cube made of Dural. One 
side of the cube has a \-in. plate-glass window 32 
in. sq. To determine if the window was sufficiently 
large, circular openings of increasing size were placed 
over the window. The aperture effect decreased to 
zero when the opening was still well within the 
window size. To determine if reflection from the 
chamber walls would affect the results obtained 
when the chamber is filled with smoke, measure- 
ments were made on a smoke by using alternately a 
black and a white rear wall. As the wall reflection 


and window size had little effect upon the measure- 
ments, the conclusion was reached that the chamber 
had been made sufficiently large. 

A blower is installed either to stir the smoke in 
order to obtain uniform distribution in the chamber 
or to exhaust the smoke after colorimetric measure- 
ments are made. 


A door serves as access to the 

















Ficourre 5 General view of apparatus, showing smoke chamly 


illuminators, optical-assembly tube, and photometers 
chamber for purposes of placing a smoke grenade in 
the chamber and firing it. To facilitate entering 
the chamber in order to clean the window. an 
additional door was constructed on the side of the 
chamber (not shown in fig. 5) 


3.2. Optical Arrangement 


The smoke is illuminated by four 150-w incan- 
descent floodlamps, mounted on a frame. The uni- 
formity of illumination on the window was checked 
by means of a foot-candle meter. 

The smoke is viewed perpendicularly through a 
lens that focuses upon a central region of the window. 
The lens is mounted at one end of a brass tube 
equipped with baffles to minimize stray light. The 
tristimulus filters are mounted in front of the 
phototubes. To achieve compactness of the viewing 
assembly, head-on-tvype phototubes were used 


3.3. Photometers 


To measure accurately the chromaticity coordi- 
nates of saturated colors, a colorimeter is required to 
evaluate low reflectance in one spectral region and 
high reflectance in another spectral region with equal 
precision. As the colors of the pyrotechnic smokes 
are highly saturated and further development will 
improve their saturation, this colorimeter requires 
highly sensitive photometers and phototubes. To 
satisfy these requirements and for expediency, 
commercially available multiplier photometers and 
phototubes were used. The Photovolt Multiplier 
Photometer, model 520—-M, known to be reasonably 
sensitive and stable, was chosen for use in the color- 
imeter. One minor adaptation of the photometers, 
required to permit the use of 5819 phototubes, in- 
stead of 1P21 phototubes for which the photometer 
was designed, was the replacement of the 10-prong 
socket with 9 dynode resistors by a 14-prong socket 
with 10 dynode resistors. To obtain direct readings 
of tristimulus values, Y, Y, and Z, on the correspond- 
ing photometers, modification of the microammeter 
circuit of each photometer was anticipated. but 

















only the Z-scale photometer required an increase in 
neter sensitivity. 


4. Results With Colorimeter 


As the source-filter-phototube combination does 
iot completely duplicate the CIE functions, best 
esults are obtained by calibrating the colorimeter 
vith standards having spectral selectivities similar 
to those of the smokes. But as no spectral data are 
ivailable for any smokes, standards of Munsell charts 
having colors similar to those of the smokes investi- 
cated were chosen and gave satisfactory results. 
Table 2 shows the Munsell notations for three differ- 


TABLE 2 


Munsell notations 


Smoke color lime 

Instrumental Visual 
, sR3 l 5R4/10 t 
” .5R 3/1 44/10 t 
4.5P3.5/14 R412 t 
D 5R 2/10 5 P2/7 t 
5Y6.5/10 SY6s t 

Ve 4 _ . 

. Y 46 4.556 f 


ent smokes reduced from data obtained with the 
colorimeter at the time of firing (f) and within 10 
min (fy) of that time, compared with notations 
obtained simultaneously by visual estimates through 
a daylight filter. The instrumental colorimetric 
data were converted to Munsell book notation, so 
that comparisons could be made with the visual 
estimates obtained by comparison with the color 
scales in a Munsell book of color. The colorimeter 
vields hues and values in substantial agreement with 
visual estimates, but vields generally higher chromas 
than visually obtained estimates. The uncertainty 
in a Visual estimate of hue and value is approximately 
one half-step. As the Munsell book does not contain 
chips with chromas as high as those of these smokes, 
the uncertainty in visual estimates of chroma is of 
the order of two steps. One possible reason for the 





discrepancy between the instrumentally and visually 
obtained chromas is that the approximation to the 
2-function is deficient in the region between 500 and 
550 mu, thereby yielding relatively low Z-values for 
yellow and red smokes. 

Table 3 shows Munsell notations of three different 
smoke colors taken on smokes fired under outdoor 
daylight (clear, sunny day) compared with notations 
obtained by reduction of instrumental data for 
duplicate smokes. The satisfactory correspondence 
between the sets of data indicates that the smokes 
were not strongly fluorescent and that the CIE 
functions had been adequately approximated 


TABLE 3 
Munsell notations 
Smoke ' 
Instrumental Visual 
Red SP3/14 SPS/14 
Yellow 7YS/12 7.5Y7/10 
Cireen 4072/3 73/3 
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Simplification of Calculations in Routine Density and 
Volumetric Determinations 


Charles T. Collett 


In density and volumetric determinations, 


the use of approximations 


plify calculations 


1. Introduction 


In the course of performing numerous density and 
volumetric determinations several methods were 
compared with regard to the calculations associated 
with them. 

Some methods involved equations that were found 
adequate and correct but too complicated to be 
desirable for routine work. For example, Gould and 
represent the density of a liquid in a 
picnometer by an equation of 26 terms and some 37 
brackets, parentheses, shillings, ete. 

Other systems were observed to lead the worker 
into dilemmas, from which he escaped by the use of 
approximations. These did not reduce the accuracy 
of his computed values, but were time-consuming and 
difficult to record in concise form. Often the steps 
of the calculations were worked out on scratch paper 
and then discarded, only the terminal values of the 
series being permanently recorded. 

The methods and equations finally selected for 
everyday use were derived from inspection of the 
component forces acting on the arms of a balance in 
equilibrium, i. e., forces resulting from the masses 
pushing downward and the buoyant effects pushing 
upward. Because these equations are based on 
considerations that are elementary to the problem, it 
would seem highly probable that they have been used 
before, and perhaps published. A search of the 
literature has not revealed them, at least not in the 
forms shown here. Certainly, they have been 
neglected in modern practice, and their simplicity 
and usefulness suggest that they should be presented 
for consideration. 


Stott glass 


2. Hydrostatic Weighing 


Consider that it is desired to find the mass, .V/,, 
of a glass object of irregular shape by the conven- 
tional substitution weighing. Assume that at all 
times a constant load, c, remains on the right pan of 
the laboratory balance, a condition that will be 
unchanged during all weighings subsequently dis- 
cussed, and the designation as right or left bemg an 
arbitrary identification without other significance. 


F. A. Gould and Verney Stott, Glazebrook’s dictionary of applied physics, 
Hil, 127 (Peter Smith, New York, N. Y.. 1990 


201 


certain methods of weighing have encouraged 


These can be avoided by minor changes in procedure that sim 
lhe use of counterpoises is shown to eliminate some buovaney corrections 


Equilibrium is established by placing on the left pan 
brass weights e and /, whose masses are MM, and M, 
and which are subject to air-buovancy effects, 2B, and 
B,. Mis the mass of c, and #2, is the air buovancy 
one. Then 


M,— B.4+-M,—B,;=M,—B 1) 


If the glass specimen, s, is then placed on the left 
pan, and equilibrium restored by removing weight /, 
M.—B,+-M,—B,=M,—B, 2) 
where \/, represents the mass of the specimen and 7? 
the air buovaney acting upon it. The 
weighing must be accomplished within a reasonably 
short time, so that the air density has not changed 
Subtracting eq (2) from eq (1), and rearranging, 


second 


M,=M,—B,-+B, 
M,- M, Vip 
py 
M,(1—**)+Vipo, 3) 


where p, is the density of the air, p,; is the density of 
the weights, and VV, is the volume of the specimen 

p, can be obtained by well-known methods,’ and p, 
is usually known. However, V, is unknown due to 
the shape of the specimen and is customarily found 
by means of hydrostatic weighing. 

A suspension is fitted below the left pan to support 
a basket immersed in a liquid, usually distilled water, 
and equilibrium is secured by the addition of weights 
q and h on the pan. Let Md, equal the mass of the 
suspension system and 2, and /’, the respective air 
and water buoy ancy effects 

M,—B,+M,—B,+ M,— B,— B’,.= M.— B.. (4) 
The next operation must be performed within a short 
period of time, that is, before the density of the air 
or the water has changed appreciably. The speci- 

N. Bauer, Technique of organic chemistry, vol. 1, part 1, Physical Methods, 
p. 273 (Interscience Publishers, Inc., New York, N. Y., 1949). N. 8. Osborne, 


E. C. McKelvey, and H. W. Bearce, Density and thermal expansion of ethy! 
ileohol and of its mixtures with water, Bul. BS @, 327 (1913) S197 (see p. 378 
































men is immersed within the basket and is subject to 
buoyancy effect B’,. Weight A is removed to restore 
equilibrium. Then 


M,—B,+M,—B,—B’,.+-M,— B’,=—M.—B,. (5) 
Subtracting eq (5) from eq (4) 
Bi = M,—M,+ B, 
M,—M,+4 M, p. (6) 
Pr 
M,—M,(1—#:). 
(1—*) 


Again an unknown, .7,, appears on the right side 
of the equation, making it impossible to use directly 
the relationship 


y 
Pw 
M,—M,(1—*) 
Pr - 
De ’ (4) 


where p, is the density of the water 

The unknown quantities may be eliminated from 
the right-hand sides of eq (3) and (7) by substitution. 
The resulting equations can be solved directly but 
contain more terms, thus: 


W,(1—**)—M,(®*)(1 pe) 
M Ps m ad Llp (8) 
| Pa 
Pr» 
Pa p. 
M,( | - my M,(1 ped 
De Pe ; (Y) 


However, if eq (9) is solved first, the value of V, may 
be used in eq (3), which is much shorter. al 

Equations (3) and (7) may be solved by the use of 
approximations, a method that has been found in 
considerable use.“ In a typical case it was found 
necessary to remove weights equal to 11.6954 g 
after the glass sample was placed on the pan. Air 
density was determined to be 0.001170 g/ml, and 
the weights were known to have a density of 8.4 
g/cm From eq (3): 


0.001170 


=; )+ V(0.001170) 


M,=11 6954 | 


11.6938-+- V,(0.001170), in grams. (10) 







For the hydrostatic weighing, equilibrium was 
maintained by removing 6.4711 g after the sample 
was placed in the immersed basket. Densities of 









* A similar equation appears in Bul. BS 9, 327 (1913) 8197, see p. 384, where 
the density of a liquid is found by immersinga sinker of known mass and volume 
‘NBS Circular 487, p. 2 (19%) 





textbook of physics, p. 175, Longmans, Green & Co., 
calculations are not mentioned 





the water and air were 0.996953 g/ml and 0.00) 
g/ml. From eq (7) 


M,—6.4711(1- oe =) 


Vs 0.996953 
M,—6.4702 
0.996953 


The first approximation is secured by using 
apparent weight of the sample in place of MW, in eq 
(11), giving 
- 11.6954—6.4702 . 

V, ae 5.2412 em’, 
0.996953 


approx. , 


This approximate value for V, is used in eq (10 


and we have 
5.2412 (0.001170) 


11.6938 + 11.6999¢ 


approx., VM, 
The second approximations are 


11.6999 —6.4702 
0.996953 


approx.s V, 5.2457 em (14 


AM, = 11.6938 + 5.2457(0.001170) = 11.6999 e 


approx 5 


These last approximate values are accepted for the 
real values. as obviously no further change will 
occur. 

A minor alteration in this procedure will simplify 
the equations and entirely eliminate the use of 
approximations. The suspended, immersed basket 
is retained, but the specimen rests on the left pan 
during the first weighing, and equilibrium is pro- 
duced by the presence of weight g. Equation (4 
then becomes 
B.. 


M.—B,+M,—B,-4 


v 


M,—B,—B,=M. (16 

For the second weighing the specimen is trans- 
ferred to the basket to be immersed as before, but 
to produce equilibrium it is necessary to add weight 


hon the pan.’ Equation (5) then is changed to 











M,—B,+M,—B,+M,—B, 
B,+M,—B,=M.—B.. 7 
Subtracting eq (17) from eq (16) 
B,—B,= M,—B, 
Vs(Pw— pa) = Mi - ats Pa 
Pr 






Pa 
M, ( oa ) 


Pw Pa 






V; 






' The transfer of sample from pan to basket is described by W. Watsor 
London, 1902), but 
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Because the densities of weight A, air, and water 
ean be secured easily, it is apparent that eq (18) 
contains no unknowns on the right-hand side, and 
the exact value of \, may be calculated at once. 
After its volume is thus determined, the object is 
weighed in air in the conventional manner, after 
vhich eq (5) can be solved 

The same glass piece previously considered was 
ilso weighed by this latter method. After the 
sample was moved from the pan to the immersed 
basket, it was necessary to place 5.2243 g on the pan, 
The densities of the air and water were respectively 
0.001170 and 0.996953 g/ml, and of the weights, 
8.4 g/cm.” 

From eq (18) 


5.2243(1—" ee =) 


— 5.2457 em’*. (19) 
0.001170 ; 


0.996953 


The sample was then weighed in air in the conven- 
tional substitution method, 11.6954 g being removed 
after the specimen was placed on the pan. Air 
density was 0.001170 g ml 

From eq (3): 


0.001170 


Vl 11.6954(1— a 


)+5.2457 (0.001170) 


(20) 
11.6999 ¢ 


3. Weighing a Picnometer by the Counter- 
poise Method 


Weighing a picnometer is a special case of de- 
termining the mass of a glass object having irregular 
shape and unknown volume, ol is therefore subject 
to the same general principles. If an attempt is 
made to use conventional substitution weighing, the 
same dilemma will of course arise, that is, neither the 
mass nor the volume is known, and neither can be 
determined without knowledge of the other. The 
volume here referred to is the volume of material 
(glass) that will displace air during the weighing, and 
does not refer to the capacity of the container, or 
so-called “internal volume.” 

Furthermore, the interior of the picnometer for 
liquids is of necessity closed to the outside except for 
limited entrances, usually capillary tubes of small 
diameter. For this reason, filling, emptying, clean- 
ing, and drying are not quickly accomplished, and 
the observer will wish to hold the number of such 
operations to a minimum. For example, a series of 
trials may be performed while the same liquid re- 
mains in the instrument, only minor adjustments of 
level being made between succeeding weighings. 

Thus the weighing of the empty picnometer and 
the weighing of the filled picnometer are often sepa- 
rated by considerable periods of time, and the air 
density, and hence the air buoyancy effects, will vary. 

However, the user of the picnometer is not 
primarily interested in the mass of the instrument; 


this value is never required except as a tare, and its 
consistent treatment as such throughout all opera- 
tions will permit the use of simple calculations with 
no approximations. By the use of a counterpoise 
the picnometer’s mass and its buoyancy in air may be 
entirely eliminated from all calculations, and in fact 
need never be known at any time.® 

The counterpoise is a second picnometer similar in 
size and shape to the first, made of glass with the 
same density, and the two instruments must be 
approximately equal in mass. A small protuberance 
is provided on each, such as the handle on the cap, 
and by grinding and polishing it is not difficult to 
adjust the masses within a milligram. However, the 
two instruments may be allowed to differ by as much 
as 0.1 g because the difference in buoyant effects 
when they are suspended from opposite balance pans 
with equilibrium maintained by weights is (A/p,—A 
Py)Pa, Where A is the difference between the masses 
of the two picnometers, and p,, ps, and p, are the 
densities of the glass, weights, and air. If A is 0.1 ¢ 
and the respective densities are about 2.5, 8.4, and 
| 0.0012 g/ml, the difference in buoyant effect will be 
negligible. 

In conventional laboratory practice the mass, 
volume of glass, and internal volume of a picnometer 
are not usually determined each time the instrument 
is used, but only at the time it is calibrated. With 
occasional checking, these values are considered true 
as long as there is no reason to suspect a change. 
In the method to be considered, ae the internal 
volume at various temperatures is determined. 

The counterpoise is placed upon the right pan 
with a constant load composed of weights ¢ and f/f, 
and the active picnometer is placed on the left pan 
with weights / and F sufficient to secure a state of 
balance, from which the equation may be written 


M,p— Bp+ Me—Be+ My— Br 


M,—B,+M,—B,4+-M,—B,;, (21) 
where ./ is the mass, and B is the air buoyancy effect 
of the object indicated by the subscript letter; 7 
and p being the active picnometer and the counter- 
poise picnometer. 

The picnometer is next filled with a sample, S, 
which will be distilled water during the calibration 
and later the particular liquid being tested for 
density. Balance is restored by removing weight F. 
Then 


M~ E+ ~ +i 


M, — B,’ + M,— Bi’+M,—B;'. (22) 

It has been observed that the two weighings may 
be performed at different times, thus a change in air 
density is to be expected. The buoyant effects in 
eq (22) will not be equal to those in eq (21) except 
by coincidence and are thus designated as B’’. 





*N. Bauer, Technique of organic chemistry, I, part 1, Physical methods, 
p. 272 (Interscience Publishers, Inc., New York, N. Y., 1949 








Subtracting eq (22) from eq (21), we have 


Assume that the small difference in mass between 
P and p in eq (21) and (22) has been compensated 
by a corresponding difference in / and ¢, so that 


Mp+M,=M,+M (24) 


It has been shown that the buovancy effects may be 
considered equal, that is, 


Also 
Bi +B =B’+B (26 


It follows that .W%, must be equal to \/,, and as 
the weights have the same density, By and B, are 
equal 


We then have 
U,—B, =M,—B;,’, (27 
which is essentially the same as eq (3 


Ms y_ yy, -M 


Ps Pr 


M,- 


where p, is the air density corresponding to eq (22 


p 
l 
) 
M, V/, . (29) 
Pa 
| — 
Ps 


As the sample in the picnometer during calibration 
is water at a known temperature, the density ps is 
read from standard tables, and the volume then is 
Ms; divided by that density. This process is repeated 
at different temperatures until a calibration curve 
has been developed for the picnometer over the 
desired range. 

When a liquid of unknown density is to be tested 
the weighing procedure is exactly the same. How- 
ever, eq (27), as shown in eq (3), can be rewritten to 
obtain 





M, M, (1—"*)+ Vso! (30) 





The picnometer having been previously calibrat: 
1’, is readily determined and 


uM. Mr (1—T )+ Visor 
Psy . 


For example, when a picnometer was calibrat: 
according to the method described, a weight 
120.9878 ¢ was on the left pan with the emp 
instrument, and 21.1467 g when it was filled w 
distilled water at a temperature of 24.288° C. T] 
densities of the weights, air, and water were det 
mined to be 8.4; 0.001170, and 0.997257 g¢/mi, respec 
tively. From eq (29) the mass of the water was 


0.001179 
8.4 
0.001170 


0.997257 


M,.=99.8411 99.9445 g. (32) 


The volume of the instrument was 


99.9445 


= 100.2194 ml. (33 
0.99% 25% 

This was used as one point on a calibration curve 

When the instrument was used to test a sample of 
heptane a weight of 116.5413 g was required on the 
left pan with the empty picnometer. It was then 
filled with the sample at a temperature of 24.992° C, 
and balance was restored by reducing the weight on 
the pan to 48.5356 g. The densities of the weights 
and the air were 8.4 and 0.001173 g/ml, respectively 
From the previously drawn calibration curve the 
volume of the instrument at that temperature was 
known to be 100.2201 ml 

Using eq (31), the density of the heptane was 


0.001173 
8.4 
ve 100.2201 


68.0057 ( an ) +(100.2201) (0.001173 


0.67964 g/ml. (34) 


WasHInecton, October 19, 1953. 
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Infrared Spectrum of Hydrogen Sulfide in the 
6,290-cm~' Region’ 





Harry C. Allen, Jr.? and Earle K. Plyler 


rhe absorption of hydrogen sulfide in the region 6,100 to 6,500 em has been measurer 


under high resolution. The rotational fine structure has been analyzed through the use of 


published energy tables for the rigid rotor 


\ classical centrifugal-distortion correction 


was applied to the rigid energy levels. It is found that there are two overlapping bands in 
this region, a normally strong, A-type band, the (nm, nm, n 1,1,1), and a normally weak, 
B-type band, which becomes strong enough to be observed by borrowing intensitv from the 


\-type band through a Coriolis interaction 
the resonating pair (2,1,0) and (0,1,2 


The B-type band is the lower component of 
The excited state inertial parameters giving 
the best fits are (1,1,1) band, A 10.398, B 


8.935, ¢ $548, pw H280.26 em 0,1,2 


band, A 10.394, B=8.918, ¢ $547, uw 6288.28 cm 


1. Introduction 


The absorption of hydrogen sulfide in the 1.6 4 
zion has been measured previously under low 
resolution and the rotational fine structure of the 
envelope analyzed [{1].° As it is now possible 
nerease the resolution 10 times, it was thought 
advisable to remeasure this region. The analysis o 
the new data confirms the main conclusions of the 
previous analysis, i. e., the absorption is due largely 
to the A-type band (1,1,1 However, the improved 
data show a second and weaker band overlapping 
the (1,1,1) band. This second band is a B-type 
band and is assigned as the lower component of the 
resonating pair (0,1,2) and (2,1,0 





2. Experimental Procedure 


The measurements of the 1.6-« band of hydrogen 
sulfide were made on an infrared grating spec- 
trometer constructed in the Bureau's Radiometry 
Section. The grating has 15,000 lines/in., and the 
rulings extend over 7% in. A lead-sulfide cell is 
used as the detecting element. Other details of the 
instrument are given in a previous paper [2]. In 
the measurement of the hydrogen-sulfide spectra at 
1.67 it is possible to use narrow slits (0.03 mm), and 
sharp lines separated by 0.15 em! could be resolved. 
Two cells were used in the measurements, one 60 em 
and the other 600 em in length. The long cell was 
so constructed that there was no loss in aperture. 
This cell transmits about 85 percent of incident 
energy. A detailed description of this cell is given in 
another publication [3]. Various pressures were used 
in the absorption cells so that both weak and strong 
lines could be observed Figure 1 shows a recorder 
tracing obtained on unruled paper for the 600-cem 


WAVE NUMBER , cm-! 


work reported herein was supported in part by the Office of Nava 
earch (under ONR C t N5ori-76, Task Order \ 
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Figures in brackets indicate literature references at the end of this paper 
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Ficure 1 Absorption bands (1, 1, 1) and (2, 0, 1 


The calculated positions of the lines are shown above the spectra. 
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for hydrogen sulfide in the region 6,100 to 6,500 cm™ 
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Ficure 2 Center of (1,1,1) band at slow scanning 8 peer 


lhe calculated positions and their intensities are placed above the spectra 


cell, with hydrogen sulfide at a pressure of 40 © 
Hg. Figure 2 represents the center of the sam 
region obtained by using a slower scanning speed 

The frequencies of the absorption peaks were ob 
tained by usinga comparison spectrum. Higher order 
of the atomic lines of mercury, krypton, argon, an: 
neon were superimposed on the absorption band 
An arrangement of shutters was used so that either 
the absorption lines or the emission line could b¢ 
recorded on the chart, and thus the shifting of posi 
tion of the atomic line by the overlapping of th 
absorption line was avoided. Several tracings wer 
measured, and the frequency given for each line j 
a result of the average of several determinations 
In these measurements the precision of most lines j 
0.1 em— or better. The experimental determina 
tions of the frequencies are listed in table 1 
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rapie 1 Com parison of observed and calculated spectra in 1.6-  regior Continued 
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TABLE 1 


Observed Calculated wave number 


Line | trite nsity wave rransition 2 ss — 
Nt number au (012 . 
cm em~* 
o4 6.0 6368. 65 s 7-4 6368. 75 0.5 
& 7 l 
5 4» 6368. 83 1.4 
10-6—9-» A368. 85 
10-6—9-s 
47 ’ 6370. 00 10-6 —9-_« 6370. 00 7 
l ) iv) 
8 7 6370. 02 3.3 
8 1.1 
O48 5, 6371. OF 
249 2 6372. 34 
2.0 6373. 39 rf t 6373. 24 1.2 
f ‘ 
0 “aya. & 
6.7 6374. 54 10 6374. 33 u 
1 10 $ 
t 6374. 61 3.0 
fi 1.0 
7 t iv ; 8 63 $1 7 
’ 8 2 
( 6375. 39 9 
il it 637 2 1.8 
1 10 € 
‘ 4.2 6376. 1 
‘ t Hf & 6376. 62 2.4 
8 8 
4 HATS. Se 
4 t 6379. ( 4, 6379. 03 
} l 6379. 46 f 
ll ; 
0 6380. 51 t f 6379. 87 5 
th 6480. 23 4.3 
l ll 6380. 69 4 
l 11 1.2 
G80. 4 Ss 
s ‘ 
s 8 63s 6381. 64 i) 
s 6381. 06 2.3 
s 8 
é 6382. 8 ; A382. &7 1.6 
( 6382. 99 8 
nA 
. ‘ 6383, 30 1.1 
" r find. 4 
” 6.( 6385. 1¢ t 4385, 23 2.2 
t 6385, 47 2.9 


3. Analysis of Experimental Results 


The previous analysis of the low-resolution data 
showed that the absorption in this region could be 
ceccounted for by an A-type band, i. e., the dipole 
moment change along the least inertial axis. The 
previously calculated spectrum was compared to the 
new data, and the strong transitions were brought 
nto agreement with the observed absorption peaks 
through the use of derivatives of the energy with 
respect to the inertial parameters [4]. Then all 
transitions of intensity sufficient to be observed were 
aleulated. The fit between the observed and calcu- 
ited spectra was very good, except that quite a 
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lated spectra in 


1.6-p region—Continued 
—— Observed Calculated wave number iii 
>, | -ntensity wave Transition oy 
number 7T o12 
cm cm 
24 4.0 6387.1 
‘ 6387.79 0.9 
265 2 6388, 03 
266 f 6388. 44 94-8 6388, 70 1.¢ 
, Rs 
11 10 6388.79 1.1 
ll 10 ‘ 
267 4.7 6389. 5 4 A389, 52 2.2 
268 1.3 6300. 47 4 6300. 49 u 
HY 2 63901. 65 
270 20 6392. 46 72-6 6302. 41 1.5 
271 0 6393. 33 
272 3 6393. Os 8 7 63903. 82 1.7 
273 5.0 6304. 07 10 , 6305. O7 1.0 
10 4 3 
274 44 6307. 55 ty— 5 6397. 73 1.1 
275 l 6398. 78 7,6 A308, 50 ° 
27 13 6399. OF 
277 4.9 6399, 88 dom 6400, 36 11 
278 2.7 M401. 2 
7v 7 (401. 80 
280 2.2 403. 45 
281 1 403. 69 
2s2 1.2 H405, 35 
283 2.0 406, OF 
M4 ; 406. 67 s ? Moe. o 1.4 
oni 1.3 6407. 67 
oN 0 408, 34 
27 4.1 (408, 77 tis- M408. 61 1.4 
28s 1.0 6409, 49 ( 6409, 50 
280 2.0 M410. 
200 20 410.60 
2v1 <2 M411. 53 
292 2.3 (412. OS 
203 1.8 (412.82 
24 1.8 M44. 21 
205 2.0 i414. 
au) 2.58 O415. 80 s 41 t 1.3 
297 1.3 417.14 
Us 2.1 Mis. 
204 2.7 6419. 15 
300 3 6421.4 
301 ~ 6423. 05 
wre 1.5 (424. 33 
WS 2.7 6425. 52 
4 ‘ 6427. 40 
WS 2.1 M428. 09 
306) 24 6429. 50 11—t 6249 10 


few of the weaker peaks were not accounted for. 
Through the use of known ground-state energies [4] 
it was possible to determine that many of these un- 
accounted for peaks resulted from a B-type band 
centered about 1 cm— lower than the main band. 
The excited state inertial parameters for this band 
were determined as described previously [1]. A clas- 
sical centrifugal distortion correction was applied to 
the excited-state energy levels of both bands [4] 

The excited vibrational state of the A-type band 
is undoubtedly the (1,1,1) [1]. The excited vibra- 
tional state involved in the B-type band is probably 
the lower component of the resonating pair (0,1,2) 
and (2,1,0). A calculation based on the quadratic 


expression for the vibrational-band system of 
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hydrogen sulfide [5] predicts this band to be at 
6,285 em™', compared to the observed 6288.30 em 
The parameters giving the best fit of the observed 


data are (1,1,1) A=10.398, B=8.935, C=—4.548, 
yp = 6289.26 cm (0.1.2) A=10.394, B=8.918, 
C= 4.547, »—6288.28 em™. 

Many of the listed peak positions cannot be 


attributed to single transitions. The overlapping 
of two or more transitions will give an observed 
peak position that will not be the true frequeney 
of any one transition. In many instances, however, 
definite assignments can be made. This is generally 
true for the extremely strong transitions because 
the observed frequency will not be greatly affected 
by weaker satellites. In the (1,1,1) band the mean 
deviation of 76 such assignments is 0.10 em™'. In 
the (0,1,2) band some 36 peaks were assigned that 
do not have a dual assignment in the (1,1,1) band 
The mean deviation in this case was 0.09 em™'. That 
so many of the individual deviations are greater than 
the mean must be attributed to the effect of finite 
slit on overlapping transitions. This is especially 
true of the (0,1,2) band, which is very weak compared 
to the main band in this region. The classical 
centrifugal-distortion correction may also be re- 
sponsible for some of the large deviations. However, 
it should be pointed out that for low J the corrections 
are small in any event, and for high J the correction 
should become classical. Another advantage of this 
method is that this correction can be evaluated 
level for level and applied to the rigid levels cal- 
culated from published energy tables. The agree- 
ment between the observed and calculated spectra 
may be seen in figure | and in table 1. The relative 
intensity scale within a given band is the same as 
that used in the previous analysis. In reporting the 
intensities, the line strengths in the (1,0,2) band 
were divided by two to give a better picture of the 
relative importance of the two bands. The P branch 
of the (1,1,1) band is roughly three times as strong 
as the P branch of the (0,1,2) band, and the R 
branch of the (1,1,1) band is approximately twice 
as strong as the # branch of the (0,1,2) band. Some 
transitions of the two bands appear to be different 
from the intensities that one would predict from 
simple theory. These intensity discrepancies no 
doubt arise from the effect of Coriolis interaction 
between the two bands. 

The new values of the inertial parameters of the 
(1, 1, 1) band agree with those determined previously 
within the estimated error. As the (0, 1, 2) band is 
somewhat weaker, and the transitions somewhat 
harder to identify, it is felt that the inertial param- 
eters for this band are not as well determined as 
those for the main band. 


4. Discussion 


In figure 1 it will be noted that the calculated 
spectrum does not account for all the observed 
shearptien lines at the high and low limits of the 
band. The calculated spectrum in these regions is 
we as eesreneiete wane, nat extendec = 
J=12, the upper limit of published ene tables 
(6). Likowtes, the other half of the rv. pair, 


| frequency region of this absorption. It 








(0, 1, 2) and (2, 1, 0), is predicted to fall in the hi 


has } ot 
been possible to locate any transitions that can |, 
definitely assigned to a band having this upper 
state; however, it would in all probability be lost jn 
the other two overlapped bands. 

In the infrared spectra of water-type molecules 
it has been observed that A-type bands are norma}; 
very much stronger than B-type bands. In {hy 
hydrogen-sulfide spectrum only two B-type bands 
(0, 1, 0) and (0, 2, 0), have been observed. It may 
be possible that other B-type bands are present i) 
the spectrum of hydrogen sulfide that are completely 
overlapped by strong A-type bands [5, 7, 8, 9]. In 
each case observed so far the overlapping bands 
have fitted the criterion for Coriolis-type interaction 
[10]. In the analysis of the rotational fine structur 
of all these bands, it has been possible to account fo 
the observed absorption by calculating a theoretica 
spectrum with effective moments of inertia. Ther 
has been no drastic perturbation of individual energy 
levels. This is in line with what might be expected 
As the inertial parameters of the two excited states 
are nearly the same, the initial separation of th 
unperturbed energy levels is very nearly constant 
Both the interaction and the energy vary nearly as 
K*; hence it is possible to fit the observed spectra 
through the use of effective moments of inertia [5 
the effect of the Coriolis interaction being absorbed 
in the inertial constants. The extent of this effect 
cannot be ascertained until the vibrational systen 
of hydrogen sulfide is better understood. The fact 
that these B-type bands are observed only wher 
overlapped by a normally strong A-type band seems 
to indicate that the intensity is interchanged quit: 
readily. 

As has been observed in other B-type bands [11 
R branch transitions between high 7 levels of adjacent 
J values are apparently stronger than would bx 
predicted by the simple rigid theory. There is stil 
no adequate explanation for this. 

The value of A=J,.—J,—TJ, lends additional sup 
port to the vibrational assignments proposed her: 
In both bands its value is large indicating that n 
equals at least 1 in each of the excited states. The 
exact value of A is of little significance because of 
the presence of the Coriolis interaction. 
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m the 
ands 7 Hermite’s n-point osculatory interpolation formula for equally spaced arguments at 
t may i intervals of h, employing the function and its derivative is very much more accurate than the 
ont it corresponding n-point Lagrangian formula and considerably more accurate than even the 
2n-point Lagrangian formula at intervals of h Also it is specially suited for interpolation 
letely in many functions (e. g., Bessel, probability) that are tabulated with their derivative. To 
|. In avoid the tremendous amount of labor in calculating the coefficients of f; and f! in the formes 
bar ds that they are usually given, Hermite’s formula is expressed as 
‘to . 
a Lot pr } x al Bhs Lea t+ Ro,(p 
ae ture i 
int to where 
retica a; a p h p 3, a, p 1 
Ther and where 
nergy 
ected r ay=K(n i } ( 2/ tay, 1 p 
states being 
of the ; 2 ni2 
stant Il p-) Il i). 
rlyas & ‘ : 
pectr: The constant k(n), which may be picked arbitrarily, is here chosen to make a, and b; integers 
la 15 The exact values of a; and 6; are given for n=2(1)11, 1 [(n 1)/2] to [m/2] so that this 
orhe formula can be applied exactly for any polynomial up to the 21st degree A schedule give 
: approximate upper bounds for the coefficients of f@% (g)h?"~A2"f(r) in R (p),. 
effect 
vsten : : 7 : 
e fact When a function f(z) and its first derivative are known at » points x,, /=1, 2, nea 
whet highly accurate interpolation formula due to Hermite is given by 
seems 
” » 
ty . ‘ ‘ , Ye ‘ ( , 
qui {((2)=2( 2 J 1—2L2"(2) (2-27) f(tdt+ DSO LE (ae) Pr — ae f'"(a R>,,(2x), 1) 
oa 
s [ll where 
jacent a " 
- , , —-_ ; 
ld be & L®(2)= IW (r—27,))/"(a r;), j=? is absent from Il’, (2) 
is still & j=1 
i and 
| sup ( " )? 
here R,,(2)= f°" (E)5 U(r4—-szy)? ((2n)!, least r,<E< greatest x 3) 
hat 7 Ua ) 
Th nn , 
en al Thus (1) is exact whenever f(x) is a polynomial of degree less than or equal to 2n—1! 
When the points z, are equally spaced at intervals of h, it is customary to alter the notation 
in z,, letting ¢ run from —[(n—1)/2] to [n/2] instead of 1 to n, where [m] denotes the largest 
Cl integer not exceeding m. Then it is convenient to choose a variable p given by r=2)+ ph and 
ren . , 7 oe or ™. . 
to let z,=2)+th. Also, f(z) =f(4o+ph) =fp=f, fled =f, and f(r) =f. Then (1), for f(x) 
h NBs considered as a function of p, is expressible as 
BS 48, ; n/2 n/2 
. S(to+ ph) : > (L®(p)}?{1—2L2"() (p—dD} fit D> LE (pp) *(p— Oks + Rap), 
i [(n—1)/2 i [(n —1)/2] 
951 
Chen (4) 
_— where now 
1952 [n/2] [n/2 
, Moe ; : ; . . - 
ae Bet L*(p)= Il (p—)) Il (i—j), (5) 
: j=—[(n—1)/2] j=—[(n—1)/2] 
and 
i9, 140 t ( {n/2] \? 
i R,,(p) ferr(ehe) MH (p—j)? A2n)!, 2 -ym SELL inp. (6) 
j j=—[(n—1)/2] ) 
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There are many advantages in the use of (1) or (4) over the ordinary Lagrangian interpola- 


tion formula given (for equal spacing) by 


2) 
f(to+ ph) LY? (p) fit R,(p), 


[(n—1)/2]} 


where 
[n/2] ’ 
R,(p)= (Hh HW (p—y)/nl, tomy SES Taya. 


j [(n—1)/2] 


[n/2] 
Thus letting (p 
j | 


j) be denoted by L“(p), and recalling that for reasonably small h 


we have approximately 
"(EA™~ A", (9) 


where A”/ is the approximate mth difference of the tabulated f(z), the remainder term for (7) 
is of the order of A"/L‘(p)/n!, whereas that for (4) is of the order of A®’"/{L(p)}?/(2n)!. 
Apart from the fact that A*"/ is usually very much smaller than A"/, the factor {1 (p)}?/(2n)! 
is equal to the square of L“"(p)/n! (where the latter is usually very small and less than unity 
so that its square is ever so much smaller) multiplied by the very small quantity »!/(n-+- 1) 

(2n). The user can appreciate the improvement by comparing the approximate upper 
bounds for the multiplier of /@"(£)/2" in R),(p) of (4) which are tabulated in the schedule at 
the end, with the approximate upper bounds for the multiplier of f° (&)A" in R,(p) of (7), 
which are tabulated in [3, p. xvi]! and from which this present schedule was calculated. Thus 
it will be apparent that (4) is a very much more accurate formula than (7). Of course, we are 
comparing (4), a confluent form of a 2n-point formula, with (7), which is only an n-point formula. 

But it is important to note than even if (4) is compared with formula (7) taken for 2n 
points at intervals of A, instead of n points, the remainder term would differ from that in (4) 
(apart from a different £ in f°"(£)) by the presence of the factor L°"(p) instead of the factor 
{L‘™ (p)}*, which has a very much smaller upper bound than the L°")(p), showing that wher- 
ever it is possible to be used, the n-point Hermite osculating interpolation formula is much 
to be preferred, regarding accuracy, to a 2n-point Lagrangian interpolation formula at 
the same interval A. This last statement becomes intuitively plausible when the osculating 
interpolation formula for » points at intervals of h is regarded as a confluent form of a 
2n-point Lagrangian formula whose 2n points lie within a range of nh so that the “average inter- 
val” between those 2” points is only half the interval of A for the usual 2n-point Lagrangian 
formula. Thus the upper bound for the remainder term of the n-point osculating formula 
would be expected to be of the order of (1/2*")th of the upper bound for the remainder term of 
the 2n-point Lagrangian formula; actual estimates show it to be even considerably smaller. 
For example, the 2-, 3-, 4-, and 5-point osculating formulas have error terms whose upper 
bounds are around 1/16th, 1/110, 1/640, and 1/3000 of the respective upper bounds of the error 
terms in the 4-, 6-, 8-, and 10-point Lagrangian formulas. 

A second important advantage in (4) is that it is suited for use with very many tables 
where the derivative of the function is tabulated alongside of the function itself. For example, 
it is useful in tables of Bessel functions of the first and second kind [4, 5, 6] which give J,(2) = 


Ji(z), Yi(z) Y;(2), and probability functions [7], and in numerous tables of more ele- 


mentary functions and their integrals, such as tables of sine, cosine, or exponential integrals 


[8, 9], where the derivative is very easy to obtain. 

However, the use of (1) or (4) in the form usually presented [1, 2] requires a considerable 
amount of computational labor which mounts considerably as the number of points zs, 
increases. It is the purpose of this present article to provide a means of using (4) with a small 
fraction of the labor involved in the direct calculation of the coefficients of f, and f. The idea 


' Figures in brackets indicate the literature references at the end of this paper, 
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behind this method goes back to a scheme first used by W. J. Tavlor for calculating Lagrangian 
interpolation coefficients for functions tabulated at real equidistant arguments [10], and which 
was generalized by the present writer for functions tabulated at nonequidistant arguments 
[11], and also for complex arguments whether in Cartesian [12] or polar form [13, 14], and 
finally even for functions that are interpolable by expressions that are not transformable into 
polynomials [15]. Recently, the writer in looking for some way to reduce the amount of work 
in using (4), observed that Taylor's idea could be extended also to the calculation of osculating 
interpolation coefficients. In place of extensive tables of the (2n—1)th degree polynomial 
coefficients of f; and f/ in (4), one requires for each separate x only some fixed quantities a, 
and b,, which are exact integers and are tabulated below. To see this, one merely expresses 


(4) as 
n/2 { A? 2L©”' (i) A\ , A? ) 
Pa h 2" (p) }* » > — ' |. ) (10 
f(ao-+ ph) L (p) a >> le “2 io—7) ys to 3) A Sie + Ran (D); 10) 
where * 
A l I’ (i )). (11) 


Now the right member of (4) or (10) without the R,,(p) gives the expression for a (2n—1)th 
degree polynomial, which, with its derivative, assumes preassigned values of f; and f{ at 2=c,, 
and moreover that polynomial is uniquely determined by the f, and f. For proof of uniqueness 
see [1, p. 85-86], where T. Fort gives a demonstration of the unique existence of a more general 
osculating formula. His proof is practically complete save for the explicit indication that 
the mode of representation of any (mn—1)th degree polynomial which is given at the bottom 
of page 85 is always possible (which is fairly obvious). Now we make use of this uniqueness 


of representation by putting f/(7)=1 into (10), so that both ff and R,,(p) are zero, f,=1, and 
we get 
L*(p) }?= Tara L sae (12) 
hi, ( A? _ 2L$ Oe) 
n—1)/2] (p i) P i 


Thus from (10) and (12), 


Ra ( ay ' b, , ah ,.) 
( ip t (p irq e's 


° Sen om ” 2) 
f(2o+ ph) ==" nat 4 + Ro,(p), (13) 
Nes ( ay 1 b, ) 
i= —(G—1)23\(Pp— 0)? © (p—?) 
where a, and 4, are given by 
a; k(n) A?, (14) 
b, k(n){ 2L' (i) A? , (15) 


and where k(n) is any suitably chosen constant of proportionality that depends only upon x3 
In the present case the k(n) was chosen as to give exact integral values for a, and 6, instead 


of rational fractional values. 
It is simplest to think of the approximation to f(z) in the concise form 


Y(af,4+-BAf,’) 
~ ¥ , 


Sa, (16) 
where 
a, =a,/(p—i)? +b (pv, (17) 
and 
B,=a,/(p—d). (18) 


? The dependence upon n of Ay, as well as of a; and 6; given below, is not indicated, so as to avoid cumbersome notation 
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In using (16), (17), and (18) with a desk calculator, it is easiest to first divide a, by p—: to 

get 8,, which is next both multiplied by Af,’ and increased by 6,. The latter, or 8,+-6,, is again 

divided by (p—?) to give a,, from which one obtains both af,and Ya, and finally S(a fh, + BAL D/Za,, 
The computation of the quantities a, and 6, was quite straightforward. Since A, 


n l ; n l 
(—1)'* ( ) (n—1)!, instead of A’, the proportional quantities ( ’ ' ,) 
iz 


i+[(n—1)/2] r+ (n 
were caleulated. Then they were multiplied by the —2Z,'"'(), which were calculated by 
differentiating the explicit polynomial expressions ,'")(p) and then setting p=7. All frae- 


n l ‘ . 
tions in 2/ 4 ) L,“'() were cleared by multiplication of these quantities, as well as 
tein er 


the (. ~~. y by some suitable integer, for each n, to vield the exact integral values 
i+[(n--1)/2)7" ' 
for a, and 6,, which are tabulated below. The a, and 6, were checked by both recomputation 
and by use in an example for every n where the answer was known exactly and where the com- 
putation by (13) (or the equivalent (16), (17), and (18)) doing the work in decimal form to 
avoid too much labor, gave agreement to 10 significant figures 

The schedule giving the approximate upper bounds for the coefficients of f°") ()h?" ~ A?" 
in the error term R,,(p), (see (4) with (6), or (13)), namely, the quantities { ‘"'(p) }*/(2n)!, 
was calculated from the approximate upper bounds for L“(p)/n! given in a schedule in [3, p 


xvi], by squaring the entries in the latter and multiplying by n!/(n+1) (2n). Thus, 
since the L(p)/n! was tabulated only approximately, in some cases to only one significant 
figure, some of the upper bounds given here for { 1‘"(p)}*/(2n)! are not at all precise. For 


example, if we take a rounded 0.01 and square it to obtain 0.0001, the true value of that square 
may be only one-fourth as large or over twice as large, depending upon whether the 0.01 was 
rounded from 0.0051 or 0.0149 Hence in the schedule below, in some cases for the larger 
number of points, only the order of magnitude of an upper bound for the { 1°" (p) }*/(2n)! is 
indicated. But more precise determination will hardly be needed there due to the extreme 


accuracy that is surely indicated even within the range of the uncertainty. 
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Schedule of approximate upper bounds for { L‘™(p 


' i" 
Figures in parentheses indicate the number of zeros between the decimal point and the first significant « 


tange of p 


Two-point 
Four-point 
Six-point 
Kight-point 
Ten-point 


0026 
5)82 
726 
10)04 
12)34 


Range of p 


Three-point 
Five point 
Seven-point 
Nine-point 


. 00021 
6)57 
S)IS 
11)60 
1306 





Kleven-point 
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Apparent Specific Volume of Polystyrene in Benzene, 
Toluene, Ethylbenzene, and 2-Butanone* 


Maurice Griffel,' Ralph S. Jessup, Joseph A. Cogliano,’ and Rita P. Park 


temperatures near 27° and 30° C of the densities 
of purified benzene, toluene, ethylbenzene, and 2-butanone (methylethylketone), and of 
solutions of a purified sample of polystyrene in these solvents. Within the precision of 
the measurements (about 1 or 2 in 10° density) the apparent specific volume of the poly 
styrene is a constant independent of the concentration in any one solvent, but depends to 
some extent upon the solvent. The lowest values, which were obtained with 2-butanone, 
are about 1 percent lower than the highest values, which were obtained with ethylbenzene 
The values are all lower than those of the virtual liquid polymer, obtained by extrapolation 
of data on specific volume above the glass transition temperature, although the data of 


Measurements have been made at 


other observers indicate that this is not true for all solvents 


In all four solvents the change 


of the apparent specific volume with temperature is less than for the liquid polymer. 
The changes in volume upon mixing of polystyrene with the solvents are shown to have 


a significant effect on such thermodynamic properties as heat, free energy, 


of mixing 


1. Introduction 


Theoretical treatments of the thermodynamics of 
polymer solutions are usually based on assumptions 
that imply that there is no volume change on mixing 
of polymers and solvents. Experimentally, however, 
it is found that there is a small but measurable 
change in volme when mixing takes place at constant 
pressure [1, 2, 3, 4, 11, 12, 15]. Because a change 
in volume will affect the experimentally determined 
values of such thermodynamic quantities as heat, 
free energy, and entropy of mixing at constant 
pressure [5], a knowledge of the magnitude of the 
change in volume is of importance in comparing 
experimental results with theoretical predictions. 
The work described in this paper was undertaken 
to provide data to aid in interpreting the results of 
measurements of heats of mixing of polymers and 
solvents. 


2. Materials 


The polystyrene used was obtained from the Dow 
Chemical Co. through the courtesy of R. F. Boyer. 
It had been prepared by thermal polymerization with- 
out catalyst at 180°C. After receipt of the material 
it was purified as follows: 150 g of polystyrene was 
dissolved in about 750 ml of benzene, and was then 
precipitated by pouring the solution slowly into 
about 3 liters of well-stirred methanol. After allow- 
ing the precipitate to settle, the liquid was decanted. 
This procedure was repeated twice. The polystyrene 
was then dissolved in 2 liters of benzene, the mixture 
was frozen at —30° C. and most of the benzene was 
removed by evacuation at temperatures not eX- 
ceeding 3° C. Evacuation was continued at a 
————— 


*This work was supported in part by the Reconstruction Finance Corporation, 
Office of Synthetic Rubber, in connection with the Government Synthetic 
Rubber Program 

Present address, The Johns Hopkins University, Baltimore, Md 
? Present address, George Washington University, Washington, D. C 
Figures in brackets indicate the literature references at the end of this paper 





and entropy 


temperature of 70° C for about 10 days. The 
temperature was then gradually increased to 120° C 
over a period of about a week, and the pressure 
over the material was maintained at less than 10~* 
mm of mercury. As a result of this treatment, the 
polystyrene was fused into a solid mass. ‘To facilitate 
subsequent handling, it was ground in a cutting 
mill, so that the final product consisted of particles 
about 3 mm or less in diameter. 

From observations of light scattering the values 
M=190,000 and B=5.2X10~* unit were ob- 
tained for weight average molecular weight and 
second virial coefficient, respectively, for the purified 
polymer.‘ 

The volatile constituents of the purified polymer 
were determined ® by heating a crushed sample to 
150° to 160° C for 2 hr in a previously evacuated 
space of known volume, measuring the resulting in- 
crease in pressure, and making a mass spectrometric 
analysis of the gas evolved. The amounts of volatile 
constituents found in this way, expressed in percent- 
age of the total weight of polymer, are as follows: 


“or 
cys 


Percentage 


Constituent by weight 


Water 0. 101 
Carbon dioxide 021 
Benzene O14 


Total volatile impurities 0. 136 


The fact that appreciable amounts of volatile impu- 
rities remained after prolonged heating and evacua- 
tion, as described above, is in accord with the find- 
ings of Grassie [10] in regard to the difficulty of re- 
moving volatile impurities from polymers. 


‘These values were obtained by M. Wales of the Polymer Structure Section 
of the Bureau 
'In the Mass Spectrometry Section of the Bureau 
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The benzene used was purified® by fractional 
freezing. The other solvents were purified by distil- 
lation. Freezing-curve measurements’ indicated a 
purity of 99.65 and 99.80 mole percent for the toluene 
and 2-butanone, respectively. The other solvents 
are believed to be of comparable purity, although no 
actual determinations of their purity were made. 
Although impurities in the solvents would affect the 
densities of both solvents and solutions, it seems very 
improbable that they would affect the apparent 
specific volume of the solute to an appreciable extent 
in view of the small differences between the values 
for apparent specific volume in different solvents. 
In table 1 the observed densities of the three hydro- 
carbon solvents are compared with “‘best’’ values of 
these properties derived from data reported in the 
literature [6] 


Tanie 1 Densities of hydrocarbon solvents 


Density 


Nita ¢ 
Solvent - 


Observed 
values 


Observed 
values 


Literature 
values 


Literature 
values 


0. Sona0 
aye 
aa 


O. 87167 
we? 
see, 


Benzene 
Toluene 
Ethylbenzene 


3. Apparatus and Method 


Densities of solvents and solutions of known com- 
position were measured by means of picnometers, 
and values of apparent specific volume of polystyrene 
were calculated from the data 

The Pyrex picnometers used (A, fig. 1) are of the 
bicapillary type described by Heller and Thompson 
{1}. The capillaries were of “precision-bore’’ quality, 
and before sealing to the picnometers they were 
calibrated by measuring at \-in. intervals the length 
of a weighed mercury thread. They were found to 
have a uniform radius of 0.494 mm. 

The completed picnometers had volumes of about 
22 cm® and were calibrated by weighing when empty 
and when filled with double-distilled water [14] at 
temperatures near 27° and 30°C. The values for the 
density of water were calculated from an equation 
derived by Tilton and Taylor [7] to represent the 
data of Chappuis [8]. Weighings were made to 0.1 
mg or better on an Ainsworth semimicrobalance, 
using newly calibrated weights. Final weighings 
were made after allowing ample time for the attain- 
ment of temperature equilibrium in the balance. 
Two picnometers were used in the measurements, 
and a third was used as a tare in weighing. All 
weights were corrected for air buoyancy on the basis 
of observations of barometric pressure, temperature 
in the balance case, and relative humidity. 


* By T. M. Mears of the Engine Fuels Section of the Bureau 
' By G. T. Furukawa of the Thermodynamics Section of the Bureau 








The positions of the two menisci relative to ref; 
ence marks on the capillaries were read to 0.01 m 
by means of a Gaertner traveling microscope. D) 
ing these readings the picnometers were immersed 
a water bath whose temperature during an expe: 
ment was maintained constant within about 0.0 
deg C. The wall of the bath through which the menis: 
were observed was of plate glass. The temperatw 
of the bath was measured by means of a strain-fr 
platinum resistance thermometer and a 
Mueller bridge. 

The solutions were made up in 60 ml bottles with 
ground-glass stoppers as follows: A clean bottle con 
taining a Pyrex-glass sealed stirring rod for use wit! 
a magnetic stirrer was weighed with its stopper 
using a similar bottle for a tare. A sample of poly 
styrene of approximately known weight was then 
introduced into the bottle through a funnel in such 
a manner as to avoid getting any polystyrene on the 
ground-glass joint, and the bottle was again weighed 
A measured amount (usually 50 ml) of solvent was 
introduced into the bottle, which was then closed 
with the ground-glass stopper. In handling the bot 
tle great care was taken to avoid getting any solution 
on the ground-glass joint. The bottle was then 
placed on the magnetic stirrer, and stirring was con 
tinued until a uniform solution was obtained. When 
mixing was complete some solvent had condensed in 
the upper part of the bottle, in the ground-glass joint, 
and the inside of the stopper. This condensed sol 
vent was removed by drawing clean dry air over these 
areas, and the stopper was replaced. The two empty 
picnometers were then weighed; following this, the 
bottle containing the solution was weighed, and the 
solution was transferred to the picnometers as de- 
scribed in the following paragraph. 

The apparatus shown in figure 1 was used for 
transferring the solution to the picnometers with a 
minimum loss of solvent by evaporation. This appa 
ratus consisted of the ‘“‘wash-bottle’’ connections 
with a male standard-taper joint that replaced the 
stopper in the bottle containing the solution, and 
with a 22-gage hypodermic needle attached as shown 
The solution was forced through the hypodermic 
needle into the picnometer by introducing helium 
into the wash bottle through the dryer (C, fig. | 
under a pressure of 2 to 10 lb/in’. Rubber bands 
were used to hold the joints together. When the pic- 
nometer was almost filled the valve controlling the 
helium supply was closed, and the pressure was re- 
leased by opening a pinchcock (not shown in fig. | 
The hypodermic needle was then disconnected 
from the filling apparatus and removed from the 
picnometer. 

The ground-glass caps (D, fig. 1) were then put 
on and the picnometers were inverted and rotated so 
as to entrap a large bubble in the bulb and thereby 
remove small bubbles adhering to the wall. The 
removal of bubbles was found to be more difficult 
with solutions of high polystyrene content because of 
the high viscosity of such solutions, and for this 
reason, the mass fraction of polystyrene in the solu 
tions was not carried above 0.1. 
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A, Pienometer; B, bottle in which 
mnections; C, tube containing drying 
ysed to force solution from bottle 

pienometer when being weighed 


i pparatus 


B 





used in filling picnometers 


solutions were made up, with “wash-botth 
went to remove water from helium gas 
to pienometer; D, ground-glass cap for 


cap containing drying agent and cotton 
et with solvent for use when picnometer was in water bath 


After the removal of air bubbles from the solution, 
the picnometers were placed in an upright position, 
any liquid entrapped in the capillaries was removed 
with a hypodermic syringe with a 27-gage needle, 


rarie 2 Er perimental values of density and apparent specific 


Mass frac 
tion poly 
tyrem 


O1VUAS 
O50127 
100008 


OOS016 
OOW0TS 
015046 
(20057 
(40027 
1OO1aS 


020076 
OS0028 
1OOLIS 


020128 
050108 
100312 


O.S71653 
STIG 
875122 
SAO TS 
SSUIO 


SES 
SHON 
SOH1253 
SA2Z152s5 
SEGUE 
SOSUNE) 
SOTOILT 


S78 CSS 


861225 
861226 
NO4N 16 
870222 


STUSN6 


of the for 


Density, g/em, 


Observed 


Cc ‘ 
culated calculated 


and the menisci were adjusted approximately to 
the proper height by means of the syringe. The 
picnometers were then placed in a holder and 
immersed in the thermostated bath to a point a little 
below the ground-glass joints, and final adjustment 
of the height of the menisci was made. Any solution 
adhering to the capillaries above the menisci was 
removed with a linen thread. While in the bath the 
drying tubes (EK, fig. 1) containing calcium chloride 
were attached to the capillaries. A wad of cotton 
in the bottom of each drying tube was wet with a 
small amount of solvent to inhibit evaporation from 
the solution in the picnometer. 

After sufficient time had elapsed for the attain- 
ment of thermal equilibrium the heights of the 
menisci relative to the reference marks on the 
capillary were observed by means of the traveling 
microscope. The picnometers were then removed 
from the bath, the drying tubes were replaced with 
the caps (D, fig. 1), and the picnometers were 
weighed after allowing sufficient time for the 
attainment of thermal equilibrium in the balance 


cause. 


4. Results 


Observed values of density of solvents and 
solutions are given in table 2, together with values of 
apparent specific volume of polystyrene calculated 


volume, and values of density calculated fram empirical equations 


ml 


at Observed, apparent 


specific volume of 
polystyrene at 


Observed 


© . ‘ 
tleulated calculated 


HKENZENE 


0 871653 x ‘ 0 
871658 
S75126 
SAO 
SAUS1I 


SHSA42 0 SORaS) 
SHAH) 5 
STINS2 

S772 

SSO144 


857471 
SOSGSS 
S5U281 
SHOZ2 
SOLLI7 
SHM77 15 
S7 5063 


ETHYLBENZENE 


0. 861230 0 
861230 
SHASI2 
870211 
s7Tvay2 


ROAST 
SSRSLS 
862136 
SOT 53S 
876753 


2-RUTANONE 


0. 708851 Ix 0 
TORSS 1 8 
SO3257 4 
SOUOTT +38 
821304 16 


7O5575 TO55ST 

TOSSS4 705587 K 

SOOU22 SOOOTY q v2 
SOWTSN SOOTII : e100 
SISICO SISI73 : vil 
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Taste 3 Values 


Benzene 
Toluene 
Ethvibenzene 
2- Butanone 


therefrom, and with values of density calculated 
from empirical equations of the form 


(1) 


where w, is the mass fraction of polystyrene in the 
solution, p is the density of the solution, and p, and 
~ are constants at a given temperature. This 
equation was derived from the definition of 
apparent specific volume on the assumption that 
this quantity is independent of composition of the 
solution and is equal to g. Values of the constants p, 
and g: were determined by the method of least 
squares from the experimental values of p and w» 
and are given in table 3. The values of g may be 
compared with the experimental values of apparent 
specific volume given in table 2. 

It may be noted that eq (1) implies that density 
can also be represented by an equation of the form 


At, 


p= Pi (2) 
where ¢ is the concentration of polymer in grams of 
polymer per cubic centimeter of solution, and where 
A (1 Pie 

Table 2 shows that the agreement between the 
observed and calculated densities within 1 or 
2 107°, except in the case of solutions in 2-butanone, 
where the maximum difference is 4 10 The re- 
sults on 2-butanone solutions are less reliable than 
the others because of the difficulty of avoiding some 
absorption of water in handling. It is believed that 
the agreement between experimental and calculated 
densities is within the accuracy of the measurements. 
This substantiates the assumption made in deriving 
eq (1) that within the accuracy of the data the ap- 
parent specific volume of polystyrene is independent 
of the composition of the solutions over the compo- 
sition range covered. This means, of course, that 
in this range the partial specific volume of the solute 
is constant and equal to the apparent specific volume, 
and that the partial specific volume of the solvent 
is constant and equal to the specific volume of the 
pure solvent. 


Is 


5. Comparison With Previous Measurements 


Values of apparent specific volume of polystyrene 
obtained in the present work are compared in figure 
2 with values for the apparent specific volume in 


of the constants pi and © 


sation 1 


em'/9 gicem' 
0. 92026 0. SORS19 
920%) S57447 
v2z2Zs81 SOS4U) 
e103 TWANG ellos 


chlorobenzene and bromobenzene by Heller a: 
Thompson [1], in benzene by Danes [11], in toluen 
by Boyer and Spencer [4], and in various solvents by 
Streeter and Boyer [15]. The values attributed to 
Boyer and Spencer were calculated from the density 
data given in reference [4]. They are somewhat 
higher than the values of partial specific volum: 
attributed to Boyer and Spencer by Spencer and 
Gilmore [12]. Streeter and Boyer reported densities 
of solutions of various concentrations, expressed in 
grams of polymer per unit volume of solution, and 
values of partial specific volume calculated from these 
data. The value of partial specific volume in each 
solvent was calculated from two experimental points 
that lie exactly on a straight line drawn to represent 
density as a function of concentration. The two 
points used for this purpose were taken in the low- 
concentration range, except where there were seeming 
discrepancies in the data. 

As stated previously, a linear relation between 
density and concentration (eq 2) implies that the 
partial specific volume is constant and equal to the 
apparent specific volume, which is related to the 
constants A and p; in eq (2) by g=(1—A)/p,. The 
use of points in the low-concentration range to de- 
termine the slope A is undesirable because, if the 
slope is constant, it can be determined more ac- 
curately by making use of the data over the entire 
concentration range. The present writers have 
found that, except for solutions in one solvent as 
noted later, the density data of Streeter and Boyer 
can be represented within 1 or 2 10~* by equations 
of the form (2), with the constants p, and A deter- 
mined by the method of least squares. There is no 
indication of a systematic trend in the deviations of 
experimental densities from such equations, so that 
within the accuracy of the measurements, the data 
of Streeter and Boyer indicate that g. is a constant 
independent of concentration for each of the solvents 
(with the exception of the one referred to previously 
For this reason, it seems preferable to calculate the 
values of g from the values of p,; and A determined 
from all the data on each solvent by the method of 
least squares, and the values shown in figure 2 were 

obtained in this way. 

The density data of Streeter and Boyer on solu- 
tions in o-dichlorobenzene depart from a straight 
line (eq 2) in a systematic manner by amounts that 
appear to be considerably in excess of the uncer- 
tainty of the data. Values of ¢g calculated from ex 
perimental values of density by means of eq (1) can 

| be represented graphically with a precision corr 
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sponding to about 1x 107* in density by a curve that 
is horizontal for values of concentration from 6 to 12 
«/100 ml of solution, but which rises at lower con- 
centrations and approaches a value of about g.—1.0 
at zero concentration. The value for o-dichloro- 
benzene attributed to Streeter and Boyer in figure 2 
was taken from the horizontal part of this curve 

Data not shown in figure 2 include a value ¢ 
0.877 at 25° C derived from the data of Streeter and 
Bover on solutions in chloroform; values of g) in 
various solvents ranging from 0.88 to 0.93 at 20° C 
with uncertainties of +0.01 to +0.04 reported by 
Signer and Gross [2]; and the values, all at 20° C, 
o — 0.939 in benzene, 0.945 in 2-butanone, 0.953 in 
evelohexane, and 0.889 in chloroform derived from 
data reported by Breitenbach and Frank [3] on the 
change in volume accompanying swelling of a cross- 
linked polystyrene containing 0.12 percent of divinyl- 
benzene. It is possible, of course, that these latter 
values are not comparable with values of g derived 
from densities of solutions of poly sty rene. 

The straight line near the top of figure 2 represents 
an extrapolation of the linear equation given by Fox 
and Flory [9] for the specific volume of polystyrene 
of infinite molecular weight above the glass transition 
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Ficure 2. Apparent specific volume of polystyrene in various 


solvents 


@. Heller and Thompson, bromobenzene; 0, Heller and Thompson, chloro- 
benzene; A, Danes, benzene; &. Boyer and Spencer, toluene; +, Streeter and 
Boyer, various solvents as follows: (1) Decalin; (2) z-chlorotriethylbenzene 

dioxane; (4) o-dichlorobenzene; (5) r-dichlorodiethylbenzene; (6) ethy! lau 
rate; (7) benzene; (8) ethylbenzene; (9) toluene; (10) ethylacetate; (11) chloro 
rm; (12) 2-butanone; (13) 2-heptanone; (14) tetralin. 

lhe circles represent results of present work as follows; ©, Benzene; @, tolu 
ne; ©, ethylbenzene; @, 2-butanone 

lhe straight line near the top of the figure is an extrapolation of Fox and Flory's 
equation for the specific volume of polystyrene of infinite molecular weight above 


the glass-transition temperature 


temperature, and may be regarded as representing 
the specific volume to be expected of bulk polysty- 
rene if it could be obtained as a liquid at room tem- 
perature [1]. It is seen that most of the values of 
apparent specific volume in solution lie below this 
line, and that all the data indicate that the coefficient 
of expansion of dissolved polystyrene is somewhat 
smaller than that of the virtual bulk liquid. This is 
true also if we include the data of Danes over the 
temperature range from 15° to 65° C, and the data 
of Boyer and Spencer from 0° to 40° C, 


6. Effect of Changes in Volume Upon 
Thermodynamic Properties 


If we take the value 0.949, em®/g calculated from 
the equation given by Fox and Flory {9] for the 
specific volume of solid polystyrene at 27° C, it will 
be seen that, according to the data given in table 3, 
the change in volume, (A1'= g).— 0.949), when 1 g of 
solid polystyrene is dissolved at 27° will have the 
following values for the various solvents: 


Solvent 


cnr’ a polystyre ne 
Benzene. 0. 030. 
Toluene . 0305 
Ethylbenzene . O28, 
2-Butanone 039, 


It is of interest to consider the effect of such 
changes in volume upon thermodynamic properties 
of polystyrene solutions.’ If polymer and solvent 
are initially at atmospheric pressure and if mixing 
is supposed to take place at constant volume, and 
then a change in volume, AV, is allowed to take place 
to bring the pressure back to atmospheric, the 
change in any thermodynamic quantity, G, due to 
the change in volume will be given by 


OG 


AG (7 


) av. 
A 


For example, we obtain from well-known thermo- 
dynamic relations 


Ak 
* Similar calculations have been made for various liquid mixtures by 
Scatchard, Wood, and Mochel [16] and by Wood and Gray [17]. 
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where @ and £8 are coefficients of expansion and 
respectively For an_ infinitely 
dilute solution in benzene at 27° C we obtain 
a 232 10-° (° C)~' from the data of table 2, and 
from International Critical Tables [13], 8 =1.00* 10 
bar~', so that at 27° C 300° K 


compressibility, 


630 jem, 


oll 
ov), 


and for the change in volume of —0.0303 ¢m* upon 
solution of 1 g of solid polystyrene in benzene 
AT =19.1 7=4.6 cal. Similarly, we obtain for the 
effect of volume change due to solution of 1 g of 
polystvrene in benzene at 27° C Al’ 2.7 cal, 
AS 0.0089 cal/° C, AF =7.2 cal. The values for 
Al’ and AH resulting from volume change upon 
mixing may be compared with a preliminary experi- 
mental value of approximately —6 cal for the heat 
of mixing of | g of polystyrene with toluene. 

In a similar manner, taking the volume of the 
liquid polymer at 27° C as 0.927, [9], the changes in 
the thermodynamic functions due to change in 
volume on mixing | g of liquid polystyrene with 
benzene to form an infinitely dilute solution are 
calculated to be A/7~=+1.3 cal, AU’ 0.7; cal, 
AS 0.0025 eal /°C, AF'= + 2.0, cal 

According to Heller and Thompson [1] only the 
change in volume upon solution of liquid poly- 
styrene, and therefore only the changes in the fune- 
tions //7, UL’, ete., caleulated from this change in 
volume have thermodynamic significance, because of 
the fact that the solid polymer is not in a state of 
thermodynamic equilibrium. It will be seen that 
even in the case of the liquid polystyrene there is a 
significant change in the thermodynamic functions as 
a result of change in volume upon mixing with solvent 

As stated previously constancy of the apparent 
specific volume implies that the partial specific 
volume of the solvent is equal to the actual specific 
volume of the pure solvent. There is then no 
change in volume due to dilution and no effect of 
change in volume upon partial thermodynamic 
quantities, such as AH,, AS8,, AP,, ete., within the 
composition range covered, 

It should be emphasized that the results reported 
in this paper are for unfractionated polystyrene of 
unknown molecular-weight distribution, which con- 
tained the amounts of impurities given previously. 
There appears to be no satisfactory method of reduc- 
ing the results to the basis of pure polystyrene. The 
results may depend to some extent upon the molec- 
ular-weight distribution in the polymer. This 
problem has been discussed by Heller and Thompson 
[1] on the basis of certain assumptions regarding the 
relation between the specific volume of the virtual 
liquid polymer (which is known to devend upon 
molecular weight [9]) and the apparent specific 
volume in solution. If it is assumed that the num- 
ber-average molecular weight of the polystyrene 
used in the present measurements is one-half of the 





O 


weight-average and that the molecular we 
distribution functions is that given by Lansing ; 
Kraemer [18], then the treatment of Heller 
Thompson [1] indicates that the measured appar 
specific volume does not differ appreciably from | 
of a polystyrene having a uniform molecular wei 
of about 100,000. 


7. Summary and Conclusions 


Values are reported for the densities of purified 
benzene, toluene, ethylbenzene, and 2-butanone, and 
of solutions of polysty rene in these solvents at 
temperatures near 27° and 30° C. Within th 
accuracy of the measurements, values of apparent 
specific volume of polysty rene derived from the adata 
are independent of composition of the solutions in thy 
range 0.5 to 10 percent of polysty rene by weight 

The results confirm the conclusion of Heller and 
Thompson [1] that the apparent specific volume of 
polystyrene depends upon the solvent. This is also 
shown very strikingly by the results of Streeter and 
Boyer [15]. 

The changes in volume upon mixing of polystyren: 
with solvent are shown to have a significant effect on 
the thermodynamic properties of the solution 
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